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ABSTRACT 


-- This thesis presents the use of a general purpose non- 
linear optimization program in the preliminary stage of ship 
design for the selection of a propeller based on methodical 
series propeller test data. The propeller series uciiizea de 
the well-known Wageningen B-Series. Three (3) “Design Cases", 
representing the thrust, power and matching approaches to 
powering problems, are formulated as FORTRAN subprogram 
analysis codes for solution by the synthesis/optimization 
program COPES/CONMIN. Designer constraints considered are: 


1) diameter limitation, 


2) cavitation limit on expanded area ratio using 
Keller's criterion icon 


3) strength requirement determined by an empirical 
relation and by a method developed by Schoenherr 
with modifications by the author. 

Objective functions considered are maximized open water 


efficiency and minimized propeller blade weight. Optimized 


solutions to specific problems previously presented by other 


authors are obtained and results are compared. 
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I. INTRODUCTION 


A. BACKGROUND 

The ship design process, in its most rudimentary form, 
has been formulated and tracked by the utilization of the 
classical design spiral (see Figure 1.1). The design follows 
a convergent helical path past each major milestone "spoke" 
until, after numerous iterative cycles, the final configura- 
tion is "centered" upon. Whether one attempts to segregate 
the principal phases of Preliminary, Advanced and Contract 
Design into separate spirals or combine these phases in 
series along the entire path to the center, it is not long 
before the designer's roughed-out sketches give way to serious 
“number crunching", specifically that of propulsion power 
estimation. 

To estimate the power required to drive the ship through 
the water at its design speed, a decision must first be made 
as to what type of propulsor (i.e., propeller, water jet, 
paddle wheel, etc.) will be used. For the average case, and 
for the discussion that follows, the marine propeller is 
chosen to be the propulsion device. Since 

a ship propeller may be regarded as a transducer 

that converts the rotational power transmitted 

through the shaft into the translational power to 

propel the ship, [Ref. 18: p. 10] 
the selection and design of this device is obviously an 
important factor in the eventual size (weight and power) of 


the ship's propulsion plant. While hydrodynamicists provide 
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a myriad of theories and techniques to generate a "custom 
built" (i.e., wake adapted) propeller for the ship under 
consideration, their expertise is usually not required in 
the early stages of preliminary design simply because the 
design has not been refined enough beyond gross estimates. 
At this stage, the designer strives to formulate what is 
possible based on previous experience. For preliminary 
power estimation, previous propeller designs (i.e., "stock" 
propellers) and results from methodical series of model 
propellers are analyzed by the designer in order to select 
the "best" available propeller under various conditions posed 
by the problem under consideration. Three examples of typi- 
cal problems encountered in preliminary ship design are: 

1) Given the ship's effective horsepower at a specific 
speed and estimates of hull performance parameters, which 
propeller, as determined by certain principal characteristics, 
will require the least amount of delivered power from the 
propulsion plant? 

2) Given the delivered power from a specific propulsion 
system in terms of torque and revolution rate at the pro- 
peller/shaft interface and estimates of hull performance 
parameters, which propeller will generate the largest effec- 
tive horsepower and speed parameters? 

3) Given a ship's effective horsepower and speed, various 
hull performance parameters, and the propulsion plant's 


delivered power characteristics, which propeller will “match” 
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these requirements at a minimum amount of weight for a speci- 
fied material? 

(Author's Note: For the sake of brevity, the three 
selection problems just cited will, henceforth, be referred 
to as "Design exe No. 1", "Design Case No. 2" and "Design 
Case No. 3", respectively.) 

For this study, the methodical propelier series method 
is viewed as the designer's choice for preliminary powering 
analysis. One of the most-widely used methodical series data 
on model propellers is the Wageningen B-Screw Series. 
Initially, the results of the series were presented as tabu- 
lations of non-dimensional thrust and torque coefficients 


(Kn and K, respectively) versus the non-dimensional advance 


Q 
ratio (J) for analytical work and as the familiar "Bp-é" and 
“Bu-6" diagrams for design purposes. As "trial & error" 
design methods performed by hand in all engineering disci- 
plines gradually transcended to numerical manipulation by 
the modern digital computer, the necessity for the adaption 
of the Series results to a format suitable for use in computer- 
aided design methods became obvious. This was accomplished 
through multiple regression analysis of the original open- 
water test data of the 120 propeller models in the Series 
and presented in the form of polynomial expressions for "Kn" 
and "Ko" (Refs. 1,2]. 

The adaptation of the Wageningen B-Screw Series polynomials 


to various types of propeller selection problems formulated 


Bessie: 


oe a ee Pe Oy ee Eee e% 
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for computer solution has been implemented recently by two 
authors. Triantafyllou [Ref. 3] and, of late, Markussen 
{Ref. 4] presented different propeller selection problems and 
proposed different schemes for computer-aided "optimized” 
solutions. In short, specific expressions for the con-~ 
straints imposed and the objective (optimality condition) 

to be maximized, expressed in terms of a number of design 
variables and parameters, were developed. Then, each system 
of equations was solved by a Newton-Raphson method to give 

a solution set of the design variables which maximized the 
objective and met all constraints. 

Rather than formulating and coding a different optimiza- 
tion scheme each time a propeller selection problem presents 
a different combination and number of design parameters, 
variables and constraints, a better approach would involve 
formulating the problem (constraints and objective function) 
once in terms of all design parameters and variables and 
utilizing a general purpose optimization scheme which can 
handle any combination and number of constraints and design 
variables. This alternative certainly allows the designer 
more flexibility in solving his problem. Moreover, it elimin- 
ates repetitive coding and debugging associated with the 
implementation of a computer-sided solution for each particu- 


lar design problem. 
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B. PROBLEM STATEMENT 

The problem, then, is that the previously cited computer- 
aided "optimized" solutions to the propeller selection problem 
are not broad enough in capability to handle variations in 
the problem formulation. The objective of this thesis is to 
apply an available general purpose optimization computer 
code to the solution of various propeller selection problems 
encountered in Preliminary Ship Design in order to enhance 


the flexibility of the selection procedure. 


C. SCOPE 

To achieve the stated objective, the general purpose 
non-linear optimization code CONMIN [Refs. 5,6] together with 
the engineering synthesis code COPES [Ref. 7] (hereafter 
referred to collectively as COPES/CONMIN) is utilized in the 
solution of the three previously cited preliminary design 
propeller selection problems. Using the Wageningen B~Screw 
Series propeller characteristics expressed in polynomial 
expressions of various design variables, three "analysis" 
codes, required by COPES/CONMIN, are developed in such a way 
that various combinations of design variables and constraints 
are used, thereby demonstrating the applicability of the 
COPES/CONMIN optimization program in the solution of propeller 


selection problems. 


D. THESIS ORGANIZATION 
The remainder of the thesis is organized in the following 


manner. 
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Chapter Il presents a short description of the optimiza- 
tion problem in general terms and a follow-on discussion of 
the COPES/CONMIN optimization program and the athena eioall 
techniques employed therein. 

Chapter III introduces definitions and concepts applica- 
ble to the propeller selection problem. A subsequent dis- 
cussion on the Wageningen B-Screw Series is followed by final 
comments on constraints imposed on the propeller selection 
problem. 

Chapter IV presents the formulation of the propeller 
selection problem as a design optimization problem which 
can be solved using COPES/CONMIN. 

Chapter V discusses the background, formulation and 
programming utilized in estimating a propeller blade's weight 
for subsequent consideration as an objective function. 

Chapter VI reviews the author's modifications to the 
propeller strength analysis developed by Schoenherr [Ref. 8] 
in the early 1960's for the American Bureau of Shipping. A 
subsequent discussion on the programming details of FORTRAN 
codes, which are utilized for the determination of adequate 
propeller blade strength, completes the chapter. 

Chapter VII reviews the formulation and programming for 
the analysis code which is used in solving propeller selection 
problems represented by Design Case No. 1. Sample solutions 
are presented and compared to those presented previously by 


other authors. 
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Chapters VIII and IX consider Design Case No. 2 and 
Design Case No. 3 selection problems, respectively, in a 
similar fashion to Chapter VII. 

Chapter X, the final chapter, presents the author's 
conclusions and recommendations. 

As a final note, all computer coding presented in this 
thesis is done in FORTRAN IV, the language used by COPES/ 
CONMIN. For the reader's convenience, Appendix A provides 
a cross-reference of the symbols presented throughout the 
thesis to appropriate FORTRAN variable names appearing in 


the author's codes. 


II. OPTIMIZATION 


A. INTRODUCTION 

The purpose of this chapter is to introduce definitions 
and concepts used in the formulation and solution of the 
general optimization problem. Then, a short discussion on 
the theory and implementation details of COPES/CONMIN is 
presented. 

For further study on the theory and methods of optimiza- 
tion, the reader is directed to the texts by Fox [Ref. 9], 


Fiacco and McCormick [{Ref. 10], and Himmelblau [Ref. 11]. 


B. DEFINITIONS 

' Before discussing the techniques of optimization and 
their application to engineering problems, some preliminary 
definitions of basic terminology should be stated. Terms 
which have relevant significance are: 

1) Parameters--The numerical quantities for which values 
are assigned to produce a design are called parameters. From 
this, it follows that a design may be specified by a vector 
D containing "p" components, each of which is associated with 


a parameter. That is: 


However, ina design process, the parameters are determined 
by some logical procedure through analysis of some kind. 
Some might take on fixed values to become "preassigned" 
parameters. Interrelationships among other parameters might 
exist so that only some of the parameters are changed when 
one design is compared to another. This consequence leads 
to the definition of "design variable". 

2) Design Variables--The parameters for which values are 
chosen in some fashion to produce a de3ien are called design 
variahles. They represent an ordered collection of components 
which is a subset of the design vector D. This subset is 
unique in that its components are “variable”, i.e., they may 
take on different values in the design process. Having 
“preassigned" or fixed some of the design's parameters and 
only allowing the remaining "desiqn variables" to change, leads 
to the conclusion that a design is now uniquely specified by 
a vector X containing "n" components (n < p), each of which 


is associated with a design variable. That is: 


oe 


: (2.2) 


ad 
" 


X 


3) Objective Function--The computable function of all or 
some of the design's preassigned parameters and/or design 
variables, with respect to which the design is to be opti- 


mized, is called the objective function. Single valued in 
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quantitative terms, the objective function's minimum or maxi- 
mum value represents the "best" obtainable or "optimized" 
design. It is expressed as F(D) to show its dependence on 
the design's parameters. But, since a design can be uniquely 
defined by X alone, then clearly F(X) suffices as an expression 
for the objective function. 

4) Constraints--Restrictions on the design which must be 
satisfied in order to produce an acceptable design are 
called constraints. A constraint may be classified as a 
“side” or a "behavior" constraint. A side constraint re~ 
stricts or bounds the range of the design for reasons other 
than direct consideration of performance. The side constraint 


on the "i"th design variable may be expressed as: 


lower x. < xupper 


Xi = i-~- L 


i= 1,...,n (2.3) 
A constraint derived from those performance or behavior 
requirements that are explicitly considered is called a be- 
havior constraint. Most often, it appears as a computable 
functional relation involving the design's parameters, both 
preassigned and variable alike. The relation may be an in- 
equality so that the "j"th of "m" inequality constraints 


can be expressed as: 


G; (D) < 0 Se a valeg ttl (2.4) 


Alternatively, the relation may be an eguality on the "k"th 


of "£" equality constraints expressed as: 


Hy, (D) = 0 Kee Lye ccyh (2.5) 


Of noteworthy importance here is that, as before, if 
some of the design's parameters are preassigned, then the 
resulting design is that defined by X which contains only 
the parameters that can be varied in the design process, 
i.e., the design variables. Therefore, constraints imposed 


upon the design may be expressed under one equation as: 


G(X) < 0 a ere 
Hy, (X) = 0 k= 1,...52 (2.6) 
ot < Xs < > ae pee eee 


A final form for constraints is that of the discrete-valued 
design variable. 

5) Feasible Design-~A design in which specified constraints 
are satisfied is called a feasible or "acceptable" design. 

6) Infeasible Design--A design in which constraints are 


violated is called an infeasible or “unacceptable" design. 


C. PROBLEM STATEMENT 
If one presupposes that a range of designs exists within 


a selected design concept, then it follows that different 
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methodologies also exist by which one may choose the param- 
eters which describe the design. One such method is optimi- 
zation where parameters are chosen in a way that the design 
will satisfy all of the limitations and restrictions imposed 
upon it and will be "best" in some sense. In view of the 
foregoing definitions, optimization is then a selection method 
applied to a design problem by which an objective function 
F(D) is minimized to produce an acceptable design which satis- 
fies a certain set of requirements called constraints. 
Formulated mathematically, the general, non-linear, 
constrained optimization problem may be stated under one 


equation as: 


Minimize: F(D) = OBJ 
Subject to: G, (D) < 0 DL pw as yi (227) 
H, (D) < 0 k = Lreweyh 
xlower 2h, a A TRPS= i 3 Saree ero | 
i — i i 


Again, as pointed out in the previous section, the design 
may be uniquely defined by just its design variables as 
specified by X when some parameters are preassigned. Thus, 
the general, non-linear, constrained optimization problem 


can now be stated under one equation as: 
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Minimize: F(X) = OBJ 


Subject to: G(X) < 0 j = 1,...,m (2.8) 
H,(X) = 0 eae a 
Hwee oi SU RPOE Gee Guan 


i i-‘i 


Solutions methods for this optimization problem are 
abundant. Those pertaining to the linear and quadratic 
optimization problems involving a few design variables are 
most often presented in graphical or analytic form, although 
numerical schemes are, by no means, a dormant form. Struc- 
tural and thermal problem solutions are most prevalent. How- 
ever, as the optimization problem becomes more complex in 
terms of non-linear relationships among an increasing number 
of design variables and of an increased number of design 
constraints, numerical or mathematical programming techniques 
dominate the solution methods. 

To limit the scope of this discussion, only the numerical 
techniques relevant to COPES/CONMIN will be considered. For 
more background on optimization techniques and applications, 
the reader is directed to a recent paper by Vanderplaats 
(Ref. 12] which presents a concise, but thorough, qualitative 
review of optimization. Although this paper deals exclusively 
with the application of design optimization to structural 
problems, it also contains a very extensive and current list 
of references on general techniques and applications of 


optimization. 
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D. COPES/CONMIN 

As previously stated in Chapter I, COPES/CONMIN is the 
collective acronym for the FORTRAN program utilizing the opti-~ 
mization code CONMIN and the synthesis code COPES. COPES 
stands for COntrol Program for Engineering Synthesis; CONMIN 
is an acronym for CONstrained function MINimization. 

1. CONMIN 

CONMIN is a FORTRAN program, in subroutine form, which 

solves the general non-linear constrained optimization prob- 


lem as stated: 


Minimize: F(X) = OBJ 

Subject to: G, (X) < 0 5 ae ee (2.9) 
xrower 2 yy < yUPPEF 5 2a o.iin 
1 = Lo 1 


Equation (2.9) applies to the entire statement. Observe 
that equation (2.9) differs from equation (2.8) in that the 
equality constraint set, given by H{X) = 0, is not specified. 
This is because the version of COPES/CONMIN used in this 
study does not consider these types of constraints. However, 
this will not pose any difficulty in solving the propeller 
selection problems previously cited. 

Again, F(X) is the objective function (OBJ). The 
vector X contains the "n" design variables (NDV). G, (X) 


are the "m" inequality behavior constraints (NCON) imposed 


lower 
i 


upper 


and Xe are the 


on the optimization problem; X 
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respective lower and upper side constraints which bound the 


“design space" over which F(X) and G(X) are defined. As 


functional relationships involving X, F(X) and G(X) may be 
implicit or explicit, but, in any event, must be continuous 
and have finite numerical values. 

When the inequality condition of equation (2.9) is 
not satisfied, i.e., G(X) > 0 for any constraint, the con- 
straint is said to be violated. If the equality condition 
is met, i.e., G(X) = 0 for any constraint, the constraint 
aS said to be active. And, finally, if the inequality condi- 
tion is satisfied, i.e., G(X) < 0, for any constraint, that 
constraint is termed inactive. Any design, defined by x, 
which satisfies the inequalities of equation (2.9) is desig- 
nated as a feasible design. Likewise, any one which vio-~- 
lates these inequalities is termed an infeasible design. 

The feasible design with the minimum objective function value, 
often referred to as the "minimum feasible design", will, 
therefore, be the optimum design. 

During the optimization process, CONMIN employs the 
Fletcher-Reeves algorithm [Ref. 13] for locally unconstrained 
problems, and Zoutendijk's method of deasible directions [Refs. 
14,15] for locally constrained problems, in a numerical pro- 
cedure which attempts to minimize the objective function, 
F(X) = OBJ, until one or more of the constraints, G; (x), 


becomes active. The numerical search procedure begins with 


an initial X vector which may or may not specify a feasible 


design. Modifications are included in CONMIN so that, if 
the initial design is infeasible, a feasible solution will 
be obtained with minimal increase in F(X). By iteratively 


updating the design vector X by the following relation: 
xiatt) 2 xla) 4 y* gla) (2,10) 


the optimization process continues by following the con~ 
straint boundaries in a direction of search § so that the 
value of F(X) decreases with each iteration q. The scalar 
a* defines the distance of travel in the direction of search 
S. The process terminates when a vector X is found such 
that no further decrease in F(X) can be made. The vector X 
is considered to be optimal and, at least, a local minimum. 

CONMIN can be used alone as a subroutine in any 
FORTRAN program where numerical optimization is desired. 
However, in order to make the optimization process more “user- 
friendly", CONMIN has been coupled to COPES in order to 
simplify its application to various types of problems. 
Further information on CONMIN can be found in previously 
cited references [5] and [6]. 

2. COPES 

COPES is a FORTRAN program that provides automated 

design and trade-off capability to the design engineer. It 


utilizes the optimizer CONMIN to provide the following six 


specific capabilities: 


1) simple analysis 

2) optimization 

3) sensitivity analysis 

4) two variable function space analysis 

5) optimum sensitivity 

6) optimization using approximation techniques 
During the execution of COPES, say for optimization, three 
principal tasks are performed: 

1) data management on the design variables and constraints 
through location assignments in a FORTRAN common block called 
GLOBCM. 

2) decision process control on the attainment of an optimal 
design vector X through multiple calls to the optimizer 
until a minimum or maximum value of OBJ is achieved and all 
G(X) are satisfied. 


+1 when 


3) evaluation of OBJ and G5 (x) at each X7 and x? 
ICALC = 2 through multiple calls to the user-provided analysis 
subprogram, SUBROUTINE ANALIZ. 

For the application under consideration in this study, only 
the optimization capability will be used. Therefore, further 
elaboration on the other capabilities is not warranted. 
Reference [7] is the user's manual for COPES/CONMIN. 
Details on the mechanics of user implementation are presented 
with subsequent illustration by example. The reader is, 
therefore, encouraged to familiarize himself with the refer- 


ence. However, at this point, it is sufficient to be aware 


of the fact that a user of COPES/CONMIN is required to: 
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1) provide a FORTRAN subroutine called ANALIZ which per- 
forms the input of preassigned parameters, the evaluation of 
the objective function and constraints during the analysis 
phase of the optimization search and the output of the 
results. 

2) provide an assembled deck of control cards required 


by COPES. 


E. CONCLUDING NOTE 

The field of optimization is both extensive and complex 
and, therefore, the foregoing presentation is, by no means, 
complete in every detail. However, it is felt that the pre- 
satus overview, in conjunction with the cited references, 
covers the necessary prerequisites that will enable the reader 
to follow the application of COPES/CONMIN to the various 


propeller selection problems in the chapters that follow. 
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III. POWERING, PERFORMANCE AND PROPELLERS 


A. INTRODUCTION 

The purpose of this chapter is to present an overview 
of the terminology and concepts that pertain to ship propul- 
sion, propeller selection and the use of model propeller 
test data. Initially, fundamental definitions used in ship 
powering problems are presented. This is followed by a 
discussion of the "classic" types of propeller design/ 
selection problems encountered by the naval architect and 
Marine and naval engineers. Propeller model testing and 
propeller performance characteristics are reviewed next. 
The chapter is completed with a discussion of the Wageningen 
B-Screw Series. 

The goal here is brevity. The reader is, therefore, 
encouraged to investigate the references cited for further 


details. 


B. DEFINITIONS 
Some fundamental terms associated with most propeller 
design/selection problems are: 
1) Effective Horsepower (P,,) -~power required to tow the 
“bare” hull (without propeller; rudder and appendage allowance 


assumed included) that generates a given resistance (Ry) at 


a given speed (V). It is determined by: 
Rp Vv 
Pp = 4-7) (3.1) 
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2) Thrust Horsepower (P,,) -~power delivered to water by 
a propeller developing a thrust force (T) and moving at a 
speed of advance (Vy) without the influence of a hull form 


ahead of it. Pr is determined by: 


Pn = 1 (3.2) 


3) Delivered Horsepower (P,)) --power delivered by shaft 
to propeller, normally specified at the outboard side of the 
stern tube. Q. is the torque delivered to the propeller; 
n, is the revolution rate of the shaft and, consequently, 
the propeller. P), is determined by: 


2a Ce Np 


p= 350 (2) 


4) Shaft Horsepower (P,)--power delivered to the inboard 
side of the stern tube having a transmission efficiency of 
Ns- Po is determined by 


Pa (3.4) 


5) Brake Horsepower (P,) --power delivered by the prime 
mover at connection flange to the power train. While P3 is 
normally associated with the prime mover's rated power at 
this connection (BHP), it can also be specified from the 


power train/propeller side as: 
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where 1p and Ne are, respectively, the bearing system and 
reduction gear transmission efficienci:s. 

6) Thrust deduction factor (l-td)--ratio of the tow 
resistance (R,) to the thrust (T) provided by the propeller. 


It is determined by 


(3.6) 


7) Wake Factor (Taylor's) (l-wt)--ratio of the speed of 


advance (VA) to the ship's speed (V). It is given by: 
Vv 
= A 
(l-wt) = + (3.7) 


8) Advance Ratio (J)--a non-dimensional value, associated 
with propeller test data presentation (see figure (3.1), 


given by the following relation: 


Vv 
V(l-wt) A 
Np D Ny D 


9) Thrust Coefficient (K,,)--a non-dimensional value asso- 
ciated with the thrust force (T) developed by a propeller 


of diameter D which is turning at a rate n, and operating 


P 
in a fluid of density p. It is defined by the following 


expression: 
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(3.9) 


10) Torque Coefficient (Ky) --a non-dimensional value 
associated with the torque (Q,) absorbed by a propeller of 


diameter D which is turning at a rate n, and operating in 


P 
a fluid of density op. It is defined by the following 


expression: 


(3.10) 


11) Open Water Efficiency (n,) --the ratio of P, to Pp 
for a propeller in open water conditions, with a uniform 
; inflow velocity field at a speed of advance Va. It is ex- 


pressed as: 


Bic. ie ect ee eee SS! ee (3.11) 


12) Hull Efficiency (ny,)--a ratio of work done on the 


ship to that done by the propeller expressed as: 


Pp Rn Vv 
a E _ _ , fi-td) 
Ny = Pe = v7, = Iowey (3.12) 


13) Relative Rotative Efficiency (np) -~the ratio of the 


actual, behind-hull efficiency to the open water efficiency. 


The value of Np does not, in general, depart from the value 


of 1.0. Most often, n, varies between 0.95 and 1.0 for twin- 


R 
screw ships and between 1.0 and 1.1 for single screw ships. 
Applicable units for the terms in the expressions above 

are: 

1) horsepower (hp) --P., Por Poe Pp and Py 

2) pounds (1bf)--R,, T 

3) feet/second (ft/sec)--V, Va 
4) foot-pounds (ft-lbf)--Q., Qp 
5) feet (ft)--D 


6) revolutions/second (rps)--n 


P 
7) revolutions/minute (rpm) --N,, = n,/60.0 
The quantities T, V,, D, Q5 and np are obtained from the 


propeller test data results. The quantities Ry and V are 
specified from the design point on the R-V curve for the hull 
under study. The quantity o is a property of the fluid in 
which the hull and propeller operate. And, finally, Ngr 1g 
and Ng are characteristics of the bearing, gear and stern 
tube systems. In preliminary design studies, nominal values, 
based on previous designs, are usually assumed unless, of 


course, these systems have been selected and actual values 


can be specified. 


C. POWERING CONCEPTS 
1. Basic Relations 
Simply stated, the fundamental powering relationship 


to be solved in ship propulsion and powering problems is: 
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= (l-td) . 
Ee ~ I=wty " "Ro Pp anno 


Utilizing the definitions just presented, equation (3.13) 


can be rewritten as: 


Rp V (1-td) 27 Qo ny 
550 ~ Tiewt) ° "RQ ° 73597 (3.14) 


Rearranging terms of equation (3.14) gives: 


R 27 Qa. n 
T  . (l-wt)V _ ; s ™p 
Tica” S50" OO ARS eg (3.15) 


And, finally, when substitutions are made, equation (3.15) 


becomes: 
TVA _ T(l-wt)v _ ams Bp ees 
550 550 Mee Mo 550 ; 
Equations (3.14), (3.15), and (3.16) provide the basis for 


different approaches to the solution of a typical powering 
problem. More background and information on the definitions 
and equations presented above may be found in Chapter VI, 
Sections 10-16 in the text by Comstock [Ref. 16] and O'Brien's 
book [Ref. 17]. 
2. Approaches to the Powering Problem 

From equation (3.16), three types of propeller selec~ 

tion problems are discernible. In the first instance, the 


propeller thrust T and the propeller's speed of advance Va 
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are taken as known quantities. The fact that T is known 
substantiates the "Thrust Approach" nomenclature given to 
this type of selection problem. In the preliminary (or, in 
some circles, conceptual) ship design phase, the specifica- 
tion of T is based upon the requirement imposed by the 
resistance of the ship (Rp) at its design speed (V) (or, the 
effective horsepower (Pp) at V) and estimates of wt and td 

in the absence of wake surveys and self-propulsion data from 
model tests. Essentially, the thrust delivered by a selected 
propeller must provide, at least, the thrust required for the 
ship hull under study. The objective in the "Thrust Approach" 
selection problem is to determine, by logical means, the 


appropriate values of Qs and n, when the open water efficiency 


P 
(n,Q) is set by the selected propeller and its performance 
characteristics. 

In the second instance, the delivered torque (Qg) and 
the propeller shaft speed (np) are taken to be known. The 
"Power Approach" nomenclature is given to propeller selection 
problems of this type because Py is known. Here, with the 
shaft and propeller speeds being equal, the torque absorbed 
by the propeller (Q,) must be, at least, equivalent to the 
delivered torque (Qg). The corresponding objective in the 
“Power Approach" selection problem is to determine, by logical 
means, the expected ship speed (V) (or, the speed of advance 
(Va)) and the associated thrust (T) that can be developed 
when the open water efficiency (n,) is, again, set by the 


selected propeller and its performance characteristics. 
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The final, and most familiar, types of propeller 
selection problem occurs when T, Vy or V, Qz and ny are all 
known. From equation (3.16), the open water efficiency (n,) 
is now established as a requirement to be met. The objective 
is, simply, to select a propeller whose open water efficiency 
(nO)e developed thrust and absorbed torque are equivalent to 
or "match" the requirements imposed. Obviously, this approach 
on the selection problem has been designated as a "matching 
problem". 

The reader is directed to the paper by Vassilopoulos 


[Ref. 2.8} for further information. 


D. PROPELLER PERFORMANCE CHARACTERISTICS 

Up until the late 1950's, much of the knowledge about 
the performance of propellers has been gained from experience 
with models. To study the relationships governing their 
behavior, a model propeller is built and run in a towing 
tank without any hull ahead of it. This is done by running 
the propeller on a long shaft projecting well ahead of a 
narrow, hydrodynamically shaped pod or “propeller boat" which 
contains the driving mechanism and recording apparatus and 
is attached to the towing carriage. The propeller advances 
into undisturbed fluid (usually water of density o and kine- 
matic viscosity v) so that the speed of advance (Vy) is known 
and the flow into the "disc" swept by the turning blades is 
uniform. For the model propeller of diameter (D) under test, 


readings of thrust (T), torque (Qg) and shaft revolutions (n,) 
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are recorded over a range of values for speed of advance (VA) 
in this "open water”™ condition. 

Using the laws of similitude, the collected data is 
reduced and scaled appropriately into the familiar functional 
relationships between the advance ratio (J) and the non- 
dimensional coefficients of propeller performance. These 
coefficients or performance characteristics, defined pre- 
viously, are: 

l) Thrust Coefficient (Kn) 

2) Torque Coefficient (Ko) 

3) Open Water Efficiency (n,) 
Figure (3.1) graphically depicts the relationship between J 
and Kn and Ko derived from test data for a propeller defined 
by a specific expanded area ratio (Ap/A,) » pitch-diameter 
ratio (P/D), number of blades (Z) and thickness-to-chord 
ratio (t/c). 

Definitions of these terms with graphical illustrations 
pertaining to various aspects of propeller geometry can be 
found in Section 15 of references [16], [17] and in van Manen's 
publication {Ref, 19]. 

More recently, highly analytical theories (lifting line, 
modified lifting line, lifting surface, etc.) for use with 
high-speed digital computers have been formulated and subse- 
quently used in "modeling", in a mathematical sense, the 


propeller and its behavior in the "wake adapted" (or, behind 


hull) condition as well as the "open water" condition. 
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Figure 3.1 Open Water Test Results--B 4-100 Series 
Propeller 


| Additional benefits derived from this approach to propeller 
performance analysis include: 

1) determination of blade section profiles along the 
propeller's blade radius (R) to achieve uniform lift and 
internal stress distributions; 

2) computation of "off-design" performance characteris- 
tics in all quadrants; 

3) subsequent determination of hull surface forces, bearing 
loads and spindle torques induced by the propeller; 

4) prediction of steady and unsteady stress distributions 
in the propeller blade using the finite element method on 
the blade of the propeller under study. 

Obviously, this approach to propeller performance analysis 
serves to: 

1) eliminate the time-consuming and expensive model 
construction and testing of propellers in tow tanks and 
cavitation tunnels; 

2) eliminate the "scaling" discrepancies which inhibit 
the reliability of design charts and model propeller data; 

3) eliminate those design charts altogether. 

As in the case with model experiments, however, the ultimate 
objective remains the same, i.e., establishing the performance 
characteristics of the propeller in terms of Kye Ko and as 
functions of J. Having these relationships enables the ship 


designer to proceed in solving the power equation (equation 


(3.13)) through any of the approaches previously discussed. 


E. THE WAGENINGEN B-SCREW SERIES 
1. Background 
The model test data of the Wageningen B-Screw Series 
have been selected for use in the powering problems to be 
solved utilizing COPES/CONMIN. The choice was driven by the 
following considerations: 

1) the Series is widely known and, despite its growing 
obsolescence, is still used in preliminary ship design 
studies. 

2) the availability of previous investigations [Refs. 3,4] 
which utilized the series, for comparative analysis of optimi- 
zation results. 

3) the applicability of the polynomial expressions for Kn 
and Ky to computer-aided analysis. 

The Series tests were conducted from 1940 through 1960 and, 

therefore, represent propeller designs (principally naval 

and merchant applications) and design philosophy of that era. 
Specifically, the Series consists of 120 model 

propellers. As is customary in methodical or systematic 

model propeller series testing, the number of blades (2), 

expanded area ratio (Ap/AQ) and pitch-diameter ratio (P/D) are 

varied systematically, while the blade outline, the profile 

of the blade's cross section along the blade radius, blade 

cross section maximum thickness (t), blade section chord 


length (c), diameter (D) and propeller hub-to-diameter ratio 


(d@/D) were kept constant for given values of Ap/Ay and Z. 


Table (I) summarizes the variations in Z and Ap/Ap for each 
set of model propellers having pitch-diameter ratios (P/D) 


of 0.4, 0.6, 0.8, 1.0, 1.2 and 1.4. 


TABLE I 


Summary of the Wageningen B-Screw Series 


Blade area ratio 


2. Series Results 

The test results of the Wageningen B-Screw Series were 
originally presented in the form of Bp-6, Bu-d and Kns Kor 
and -J diagrams. As stated in Chapter I, multiple regression 
analysis was performaed (again, [Refs. 1,2]) on the results 
to produce the polynomial expressions for Kn and Ko: The 
open water efficiency (ng) as a function of J, follows from 
equation (3.11). The correction for "scale effects" was 
achieved by using Lerb‘s method of equivalent profiles 
(Ref. 20]. Although Triantafyllou's thesis [Ref. 21] sug- 


gests an improved method for scale correction, the results 


of Reference [2] will be used in this study. 
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In propeller selection problems which use this 


Series, values for Ky and Ko are defined as: 


an 
\ 


£,(3,P/D,Ap/Ap,2,t*/e sep) 


(3.17) 


x 
" 


£5(5,P/D,Ap/Agr2,t*/C 6p) 


To compute Kn and K the following equations are used: 


Q’ 
Ky = K! + AK 

(3.18) 
Ko = Ko + AKo 


The polynomial expressions found in Tables (5) and (6) of 


Reference [2] are then used to evaluate to components Kae 


K;, AK 


Q and AK,. 


T Q 
Table (5) in Reference [2] lists the coefficients used 


in the polynomial expressions for K,, and K! at an equivalent 


T Q 
Reynolds number (Rn 75R) of 2 x 10°, It is defined as: 
/ 2 2 
ie 7 Cc 75R (v,) + (0.751n,D ) a 
.75R v i 


where: 
v = kinematic viscosity of the fluid (£t2/sec) ; 


= blade section chord length at 3/4 propeller 
radius (.75R) in feet (ft). 
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To account for "other effects", coefficients AKn and AK, are 
introduced. Table (6) in Reference [2] lists the coeffi- 
cients used in the polynomial expressions for these coeffi- 
cients. "Other effects" include the operation of the 
propeller at an equivalent Reynolds Number different from 

2 x 10°, Also, variations in other parameters which define 
the propeller's geometry, specifically t/c values different 
from the ones fixed by the Wageningen propellers, are taken 
into account by corrections to the equivalent Reynolds number. 
By keeping the blade section's chord length (c) at the value 
of the Wageningen propeller, a change in a blade section's 
maximum thickness from the standard one defined by the Series 


(t ) to one preferred in the selection (t* ) produces a 


new equivalent Reynolds number (Rn* 5.2) given by: 


" 1+2(t/c 75R) 
Rn 35p = ©XP [4-052 + Tr2(t*7e 7.0) (2n Rn 75p74- 6052) | 
(3.20) 
where: 
Rn” oe, = the new equivalent Reynolds number; 
(t*/c) 75R = new equivalent t/c at 3/4 propeller 
s radius; 
Rn = the Reynolds number computed by equation 
.75R 
(3.19); 
(t/c) yep = standard equivalent t/c at 3/4 propeller 


radius for Wageningen propellers. 
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For the Wageningen B-Screw Series, the standard equivalent 


t/e is given by: 


(0.0185 - 0.001252)2 


t/¢ a5R = 2.073 A,/A, (3.22) 


Further details on blade section geometry will be addressed 
in Chapter V. Reference [2] contains background and other 
information on the equations above. 
3. Limitations on Series Data 
In utilizing the Wageningen B-Screw Series in any 
propeller selection problem, the following restrictions 
apply to the Series data: 
1) Number of Propeller Blades (2)--The Series considers 
only propellers with numbers of blades as shown in Table 


(I). Therefore, 
Z = 2, 3, 4, 5, 6 or 7 (3.22) 


However, the two bladed propeller, i.e., Z = 2, is not very 
common in conventional merchant and naval ship designs and, 
therefore, is not included in this study. 

2) Equivalent Reynolds Number--The Series data, as pub- 


lished, is valid only in the range of equivalent Reynolds 


numbers given by: 


If the equivalent t/c is varied from the standard equivalent 


value (t/c) 75R! then the new equivalent Reynolds number 
(Rn* sen), which results from this variation, must lie within 


the same limits. That is, 
2x10° < Rn < 2x10? (3,24) 


These limits are appropriate for full-size propellers. For 
example, given the following: 

a) wt = .22 
= 20 (knots) 


P 


Vv 
c) N, = 104 (rpms) 
D = 25 (ft) 


e) C 25pR = 4.0 (ft) 
£) v = 1.228510 (£t*/sec) 


the value for Rn is equal to 3.619 x10’, 


~75R 
3) Pitch-diameter Ratio (P/D)--The series data considers 


only pitch diameter ratios in the range given by: 
0.4 < P/D < 1.4 (3.25) 
4) Advance Ratio (J)--An inspection of the Series results 
in graphical format shows that J varies over a range given 


by: 


O: too) 156 (3.26) 


49 


5) Expanded Area Ratio (A,/A,) --Using Table (I), A,/Ap 
varies over certain ranges depending on Z. This is stated 


as; 


0.35 < A,/A, < 0.8 Z= 3 (3.27) 
0.40: < Aa. < 1,0 Zn 4 (3.28) 
0.45 < AA, < 1.05 255 (3.29) 
0.50 < A,/A, < 0.8 Z = 6 (3.30) 
0.55 < AJA, < 0.85 227 (3:91) 


6) Hub diameter-to-Propeller Diameter Ratio (d/D)-~From 
Table 37, Section 17 of Reference [16], the Series data 


requires that: 


qd/D = 0.18 Z 


{l 
Le) 
~ 
~J 


(3.32) 


ad/D = 0.167 Z 


WW 
ry 
- 
nn 
. 
a 


(3.33) 


F. SUMMARY 
From the preceding discussions, the following observa- 
tions can be made: 
1) the "Design Cases", defined in Chapter I, are examples 


of the powering equation solution approaches. That is, Design 
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Case No. 1 constitutes a "Thrust Approach" problem; Design 
Case No. 2, a "Power Approach" one; Design Case No. 3, a 
"Matching" problem. 

2) equations (3.17) and (3.8) imply that an optimization 


solution to the "Design Cases" will involve P V, wt, D, 


B’ 
* A : : 
P/D, Ap/Ao» (t/c) 5p, Q, and n, as possible design variables. 
3) when viewed from the concepts on optimization presented 
in Chapter II, equations (3.25) through (3.31) constitute 
side constraints to an optimized solution of a propeller 


selection problem which uses the Wageningen B-Screw Series. 


Having noted these points, the propeller selection problem 


can now be formulated as a general, non-linear, constrained 


optimization problem. 


IV. PROPELLER SELECTION--AN OPTIMIZATION PROBLEM 


A. INTRODUCTION 

The purpose of this chapter is to present the formula- 
tion of the propeller selection problem as an optimization 
problem that can be solved using COPES/CONMIN. Three usual 
restrictions considered by the designer in any propeller 
selection analysis are stated as constraints. Then, the 
components of D and X are assembled based on requirements 
from previously cited relationships. The restrictions 
considered by the designer and the limitations imposed by 
the use of the Wageningen B-Screw Series are presented in 
inequality constraint format. A formal statement of the 
propelier selection problem as an optimization problem is 
followed by a review of the GLOBCM common block format and 
the basic subprograms used in all three versions of SUBROUTINE 
ANALIZ that pertain to each Design Case. 

The FORTRAN subprogram listings are found in Appendix 
B. Comment cards have been used extensively in the coding 


development to assist the reader. 


B. DESIGNER'S CONSIDERATIONS 


1. Propeller Size 


The first restriction on the selection of any 


propeller is size. That is, the propeller race in the stern 


of the hull under consideration will only accommodate a 


propeller of some given maximum diameter (Dy im): 


straint on a selected propeller of diameter D, this may be 


written as: 


lim 
or, alternatively, as: 
G(X) = po -1 < 0 
De. - 
lim 
2. Cavitation 


As a con-~ ‘ 


(4.1) 


(4.2) 


Another item of importance in propeller selection is 


the cavitation phenomenon. 


When a propeller of given diameter 


D and expanded area ratio Ap/Ag is operating to produce a 


thrust T, the formation and subsequent collapse of water vapor 


bubbles on the blade surface, i.e., cavitation, is likely to 


occur if the localized surface pressures, usually 
"back" side of the blade, drop below the pressure 


the fluid would boil (p 


watvap) in the surrounding 


Avoidance of cavitation can be reasonably assured 


a propeller having certain geometric characteristics. 


on the 

at which 
environment. 
by selecting 


A good 


empirical relationship that establishes these characteristics 


for propellers typified by the Wageningen B-Screw 


the Keller Cavitation criterion [Ref. 2: p. 259]. 


fies the minimum required expanded area ratio (A, /A 


avoid cavitation and is given by: 
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Series is 
It speci- 


(0) min) 2 


where: 


D 


Patm 


Pwatvap 


As a constraint on the propeller selection, this requirement 


is writen as: 


or 


3. Strength 


The final designer's consideration (for this 


included in the selection of a propeller, 


+ epee 


(1.3 + 0.32) T 


(Patm + Macgh., ~P 


: +b 
watvap? p2 


number of blades; 
developed thrust (libf); 
propeller diameter (ft); 


atmospheric pressure (psia); 


fluid vaporization pressure (psia) 


fluid density (1b£ sec*/ft‘) 


32.174 (£t/sec*) 
depth to shaft centerline (ft) 


constant: 0.1 for Z = 2, 
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0.2 for 2 


(4.3) 


=l. 


(4.4) 


(4.5) 


study), 


is that of strength. 


Given the propeller's material (promat), selected from Table 
(II), and the loadings (T and Q,) imposed, it is important to 
ensure that the blade's cross sections have proper dimensions 
(in an ideal sense, maximum blade section thickness (t*) and 
chord length (c)) to ensure adequate strength. Since the 

use of the B-Screw Series requires that the chord length 

(c) vary as a prescribed function of D, Z and AL/A , as given 
in Table 1 of Reference [2], the adequacy for strength can 

be determined by an appropriately selected value for blade 
section maximum thickness-to-chord ratio (t*/c) alone. So, 


if CF in is the established minimum blade section maximum 
thickness, then the strength requirement follows from the 


constraint given by: 


t* 


© 5 / ee ae he) < t*/e (4.6) 


min 

The fact that blade section maximum thickness for the B-Screw 
Series varies linearly with the propeller radius (R) allows 
the strength constraint (equation (4.6)) to be evaluated at 
one section along the radius. This point is chosen to be at 


the 3/4 radius (.75R). Therefore, equation (4.6) becomes: 


(ete) pt. ENE) Gas (4.7) 


or 


> (t*/¢) 75R min 
Gy (x) = aes 6S.) era -l1 < 0 (4.8) 


Reference [2] suggests the following empirical relation for . 
the minimum required equivalent blade section maximum 


thickness-to-chord ratio (t*/c) 75R ea) 


(2375-1125P/D) P, 


Z ¢0.0028+0.21 


2 w2 
4.123N,D°(s_ +? Np ) 
* _ 12.788 
(€*/¢) 35R min = es. Aon 
° E’ "Oo 
(4.9) 
where: 
D = propeller diameter (ft); 
Py = delivered power (hp); 
N, = propeller revolution rate (rpm); 
S, = propeller material allowable stress (psi); 
P/D = pitch-diameter ratio. 


However, in Chapter VI of this thesis, an algorithm which 
employs the Schoenherr formulation [Ref. 8] with some modi- 


fications, is presented as an alternative to equation (4.9). 


C. THE DESIGN VECTOR 

In view of the preceding presentations on optimization 
and powering, the design vector DB can be assembled for the 
general propeller selection problem utilizing the B-Screw 
Series. This vector is composed of preassigned parameters 


relating to environmental conditions, hull characteristics 
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and the propeller which are required for various equations 
and the design variables. 
1. Parameters 
a. Environmental 
These parameters pertain primarily to the fluid 
conditions in which the propeller operates and to the atmos- 
phere. Required for various calculations, they are: 
1) fluid temperature (°F) --Temp 
2) fluid density (lbf sec?/ft4)-~ 
3) fluid viscosity (£t*/sec) --v 
4) fluid vaporization pressure (Psia)--Pyatvap 
5) atmospheric pressure (psia)--P iin 
b. Hull Characteristics 
These parameters pertain to certain details pre- 
scribed for the hull under study in the powering analysis. 
They are: 
1) wake fraction--wt 
2) thrust deduction--td 
3) relative rotative efficiency “"Np 
4) number of propellers-~noscrw 
5) shaft centerline depth (£t)--ho, 
6) propeller diameter limit (£t)--Dy iy 
c. Propeller 
These parameters are specified in view of their 
discrete-valued nature. They are: 


1) number of blades--Z 


2) material~-promat 
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Table (II) lists materials and properties considered in this 
study. These values are taken from Table (35), Section 15 


of Reference [16]. 


TABLE II 


Material Identifier Reference 


promat Material Allowable Stress--Sc Density--wd 
(psi) (Ib£/in?) 
1 Cast Iron 3600--3950 .260 
2 Cast Steel 5915--6265 289 
3 Type 2 Bronze 7200-7585 305 
4 Type 4 Ni-Al 8910--9430 278 
Bronze 
5 Stainless Steel 5400--5500 283 


2. Design Variables 
In view of equations (3.13) through (3.17), the 


design variables common to all selection approaches are: 


1) P 


E 
2) V 
3) 
4) P/D 
5) A,/A, 
6) (t*/c) J6p 
7) Np 
8) Qs 


The vectors D and X are shown schematically in figure (4.1). 
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Figure 4.1 
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Design Vectors D and X 
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D. CONSTRAINTS 

Besides the constraints imposed by equations (4.2), 
(4.5) and (4.8), equations (3.25) through (3.31) are rearranged 
to the format of constraints in equation (2.9). They are 
listed as follows: 


1) Equivalent Reynolds Number--Equation (3.25) becomes: 


l1< —- = 1000 (4.10) 
TO 


Two constraints are derived: 


Rn* 
G(x) = 1- -75R 0 (4.11) 
2x10 
* 
Rn 
G,(%) = -75R _ 1900 < 0 (4.12) 
2x10 


2) Expanded Area Ratio--Equations (3.27) through (3.31) 


become: 


(Ag/Ay) (2) < Ap/Ay < (Ag/Ap) (2) 


lower 


Two constraints are derived: 


Gy(X) = (Ag/AG) rower (2) - Ag/Ay < 0 (4.14) 
Gg(X) = Ap/Ay - (AB/AQ) upper (2) < 9 (4.15) 
60 
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3) Advance Ratio--Equation (3.22) becomes: 


9 < resi (4.16) 


G(X) = ~% < 0 (4.17) | 
Go) = ee T2010 (4.18) ! 
2 ~ Te = : | 


4) Equivalent Blade Section Maximum Thickness~to-Chord 
Ratio--Using equation (3.19), boundaries on the range of (t*/c) 7I5R 


are defined by: 


$(t/e) aap < (t*/C) gen < 4(t/C) g5p (4.19) 


Two constraints are derived: 


A 
[>] 


(4.20) 


x 1 
Sh BOOS) age. SEAZS) eR 


Gg (X) = (t*/c) o5p = 4(t/c) 6p < 0 (4.21) 


E. OBJECTIVE FUNCTIONS 

Upon consideration of equation (3.13), Design Case No. 1 
and Design Case No. 2 require that the open water efficiency 
(n,). given by equation (3.11), be maximized. In terminology 


related to optimization, this is stated as: 
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OBJ = - 7 (4.22) 


Design Case No. 3, the "matching" problem, requires that the 


blade weight (bldwt) be minimized. This is stated as: 
OBJ, = bldwt (4.23) 


F. PROPELLER SELECTION OPTIMIZATION PROBLEM STATEMENT 
As a general, non-linear constrained optimization problem 
to be solved by COPES/CONMIN, the propeller selection problem 


for all Design Cases may be stated as one equation given by: 


Minimize: F(X) = OBJ, » or OBJ, (4.24) 
, 
Subject to: G, (X) < 0 ee eee 
COTE ee RIPPERS dee 0098 
The constraint G2 (X) and the values for ylower ana xUPper 


i i 


will be specified according to each Design Case. 
G. CODING FUNDAMENTALS 
1. GLOBCM Common Block 
The GLOBCM common block, required by COPES/CONMIN, 
is now assembled. Table (III) specifies the assignment loca- 
tions for the FORTRAN variables which define objective 
functions, design variables and constraints. 
2. SUBROUTINE ANALIZ 
While each Design Case uses a different approach, 
all analyses are very similar. Therefore, each SUBROUTINE 
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TABLE III 


Global Common (GLOBCM) Catalog 


Global FORTRAN DEFINITION 

Location Name 

1 ETAO No 

2 WEIGHT bldwt 

3 AEDVAO A,/Ap 

4 DIA D 

5 N N, = 60 - ny 

6 PE Pe 

7 PDIVD P/D 

8 Qs Qs, 

9 TC75R (t*/c) 5p 

10 Vv V (ft/sec) 

11 RJCHL G, (X)--eqn (4.17) 

12 RJCNU G, (X) --eqn (4.18) 

13 R75RCL G(X) --egn (4.11) 

14 R75RCU G,(X)--eqn (4.12) 

15 AEAOCL G,(X)--eqn (4.14) 

16 AEAOCU G_(X)--eqn (4.15) 

17 TC75CL G. (X) -~eqn (4.20) 

18 TC75CU G, (X) --eqn (4.21) 

19 POWBAL Gj 2 (X) --eqn (7.10) or (8.11) 
or (9.4) 

20 DIACNU Gg (X) --eqn (4.2) or (9.6) 

21 AEAOCV Gy 9 (X)--eqn (4.5) 

22 TCSTRS G, , (X) --eqn (4.8) 


J 


ANALIZ shares a common structure and other common subroutines 
which perform calculations required in all cases. Appendices 
C, F and I contain, respectively, the source listings of 
SUBROUTINE ANALIZ for Design Case No. 1, Design Case No. 2 
and Design Case No. 3. 

a. Structure 

The structure common to all cases follows 

accordingly: 

1) all initialization of environmental, hull and propeller 
parameters is accomplished in the input section (ICALC = l). 

2) evaluation of Ky and Kos all constraints and appropriate 
objective functions As or bldwt) are accomplished in the 
execution section (ICALC = 2). 

3) output of results for each optimization problem is 
accomplished in the output section (ICALC = 3). 

b. Basic Subprograms 
The following FORTRAN subprograms are used in 
all three SUBROUTINE ANALIZ codes: 

1) SUBROUTINE CH75RA--calculates the equivalent blade 
section chord length (¢ o5p) for the propeller using Table 1 
(Ref. 2, p. 252]. 

2) SUBROUTINE REY75R--calculates the equivalent Reynolds 
number (Rn* 5p) using equations (3.19) and (3.20). 

3) SUBROUTINE COEFSA--calculates the thrust and torque 
coefficients (K, and Kp) through sequential calls to SUBROU- 


TINE CALCKT and SUBROUTINE CALCKQ. The polynomial expressions 
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(Tables (5) and (6), (Ref. 2]) for these coefficients are 
contained in SUBROUTINE CALCKT and SUBROUTINE CALCKQ respec- 
tively. 

4) SUBROUTINE OPWEFF--calculates the open water efficiency 
(n,) using equation (3.11). 

5) SUBROUTINE JCNA--calculates the constraints on the 
advance ratio (J) given by equations (4.17) and (4.18). 

6) SUBROUTINE REYCNA--calculates the equivalent Reynolds 
number constraints given by equations (4.11) and (4.12). 

7) SUBROUTINE EXTCCN--calculates the constraints on ex- 
panded area ratio (Ag/A) and equivalent blade section maxi-~ 
mum thickness-to-chord ratio (t*/c a6) given by equations 
(4.14), (4.15), (4.20) and (4.21). 

8) SUBROUTINE DICNUA--calculates the constraint on the 
propeller diameter (D) given by equation (4.2) using the 
hull parameter on maximum diameter (Dyin): 

9) SUBROUTINE CAVCNA~-calculates the constraint for 
cavitation given by equation (4.5) using equation (4.3). 

10) SUBROUTINE STRCNA--calculates the constraint for 


strength given by equation (4.8) using equation (4.9). 


H. SUMMARY 

The propeller selection problem has now been formulated 
as a constrained optimization problem which can be solved 
by COPES/CONMIN. Two items remain for discussion before 
proceeding to specify the final details pertaining to each 
SUBROUTINE ANALIZ code and to present numerical examples. 


These items are: 


1) the theory and coding relating to the computation of 
the propeller's blade weight for the evaluation of the objec- 
tive function in Design Case No. 3 (OBJ,). 

2) the theory and coding relating to the computation of 
the minimum required equivalent blade section maximum thickness- 


to~-chord ratio (t*/c ) for use in the alternative 


~75 min 
evaluation of the strength constraint given by equation 


(4.8). 
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V. PROPELLER BLADE WEIGHT--AN OBJECTIVE FUNCTION 


A. INTRODUCTION 

In this chapter, the method for the computation of the 
propeller's blade weight (bldwt), the objective function 
OBJ, is examined. First, a brief overview on the steps in 
the computational procedure is presented. The FORTRAN 
subprogram SUBROUTINE WGTCAL developed from the algorithm is 
then described. Again, Appendix B contains all subprogram 


listings. 


B. THEORY AND PROCEDURE 

Given the material of a propeller, the calculation of 
the weight of one blade involves nothing more than a volume 
calculation, a relatively routine task performed by most naval 
architects/marine engineers. Analagous to the determination 
of the underwater volume of a ship's hull, the calculation 
is an integration of blade section profiles' cross-sectional 
areas over the propeller radius (R). 

1. Limits of Integration 

Figure (5.1) depicts a side elevation view of a 

blade and hub, parallel to the propeller shaft axis. The 
cross-hatched area indicates the trace of the volume to be 
calculated. In view of equations (3.32) and (3.33), limits 
of integration are from r = .167R to r = R for Z2 = 4,5,6 and 


r = .18R to r = R for Z = 3,7. For convenience, a non- 


dimensional variable "x" will be defined as: 


a 
Be) AD 
™ a 
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fea) - 
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Figure 5.1 Propeller Blade & Hub--Side View 


(5.1) 


ah 


Limits are now expressed as x = .167 or .18 to x = 1.0. 
Quite obviously, R= (D/2.0). 

2. Blade Section Profile 

Figure (5.2) depicts an expanded cylindrical blade 

section in profile view at a given r or x. For the Wageningen 
B-Screw Series, the profile is defined, geometrically, by 
a succession of vertical ordinates which specify points along 
the blade section's profile on the "face" (Y¢) and on the 
back (y,) with respect to the pitch reference line. At any 
‘yr = XR, vercical ordinates for "aft" (P < 0) and "fwd" 


(P > 0) portions of the blade section are determined by: 


Ye, “= Vy Ce" = tee) P <0 (5.2) 

Yoa = (Vy tVa) (t*~ te) + ti, P< 0 (5.3) 
and 

Yer = Vz,(t* - th.) P > 0 (5.4) 

Yor = (V,+VQ)(t*- tt) + th, PR > 0 ($.5) 
where: 

Vi1V5 = tabulated values depending on x and P 


(see Tables (2) and (3), Ref. [2]); 
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tee 


it 


blade section leading edge thickness (ft); 
tt. = blade section trailing edge thickness (ft). 
Units for Yea' Ypa’ Yer and Yp¢e are feet (ft). For this 


study, a reasonable assumption is made in that: 
af 1 
t = t = (rm) t* (5.6) 
3. Blade Section Cross-Sectional Area 


The cross-sectional area at each x = r/R (A(x)) is 


determined by: 


; 
De Ait Ae; (5.7) 
i=l eo 1= fi 
¢ 
where 
h + 
3 aitl ai 
LO iL 0 ei, rs (5.8) 
Aoi Ypai ~ Yeti (5.9) 
Bt,  Ypaiay * Yeaiet (5.10) 


Expressions for Agy he; and he follow in similar fashion. 


itl 
Values for Yeai and Ypai ore determined at 9 points 


along the "aft" portion of a given blade section's chord (c): 


values of Vegi and Ypei ore determined at 10 points along the 


"fwd" portion. The values for AP Ai and APs; are fractional 
values of the blade section's chord length (c) at radius 
r= xR as determined from Tables (2) and (3) in Raterencs 
{2]. The units for A(x) are square feet (£t2) . Units for 


hie La Nei Nesay? Cc, APLi and AP ey are feet (ft). 


4. Volume Integration 
The blade volume (bldvol) is finally determined by 


using Simpson's Rule for integration of A(x) along the non- 
dimensional radius x using appropriate limits. 
5. Blade Weight 
Once the blade volume (bldvol) is calculated, the 


weight (bldwt) is determined by: 
bldwt = bldvol - wd - 1728 (5.11) 


where: 


bldvol volume of one blade (£t3); 


wd = material weight density (1b£/in?) . 


Weight Density (wd) depends on blade material (promat). 


Table (II) lists appropriate values. 


C. CODING 
SUBROUTINE WGTCAL is the main subprogram for the blade 
weight calculation. It, in turn, calls the following FORTRAN 


subprogram for various calculations: 
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1) SUBROUTINE TDIST--generates, at specified radius 
values, a distribution of blade section maximum thicknesses 
(t*). 

2) SUBROUTINE BLDPRP--generates, at specified radius values 
(i.e., vr = .167R or .18R, .2R, .3R, .4R, ..., .9R, 1.0R), 
various “blade section properties", one of which is a blade 
section's cross-sectional area given by equation (5.7). 
Other properties which are determined (for later use in 
direct stress computations) include blade section chord 
lengths and centroids and "critical point” locations as 
defined in Chapter VI. 

3) SUBROUTINE BLDVOL --performs a Simpson's Rule integra- 
tion for the propeller blade volume (bldvol) using blade 
section cross-sectional areas generated in SUBROUTINE BLDPRP. 
The blade weight (bldwt) is computed as a final step in the 
Main subprogram SUBROUTINE WGTCAL. 

Examination of the codes in Appendix B reveals extensive 
use of common blocks for passing data from one subprogram to 
another. Comment cards provide a full definition of all 
common blocks as well as a description of the task being 


performed at various points in a given subprogram. 


D. SUMMARY 
The coding developed for this study is, admittedly, not 
very compact and efficient. However, the intention has been 


to write all codes with sufficient documentation in order to 


facilitate the reader's understanding of the algorithms employed 


as well as to make the author's debugging work easier. 
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VI. THICKNESS-TO-CHORD RATIO--A DESIGN CONSTRAINT 


A. INTRODUCTION 
The purpose of this chapter is to examine the development 

of an algorithm that will be used to determine the minimum 
required equivalent blade section maximum thickness-to-chord 
ratio used in equation (4.8). The formulation is based upon 
the method developed by Dr. Karl E. Schoenherr [Ref. 8] in 
1963. After a review of the past and present methods employed 
in propeller strength analysis is conducted, a description of 
the Schoenherr model and a list of the assumptions used with 
that model is presented. Then, a brief restatement of his 
model's equations which are used in the algorithm is followed 
by a derivation of the author's modifications to the Schoenherr 
method. The chapter is completed by conducting a review of 

the theory and coding employed by the algorithm. 

The principal reference which is cited throughout this 
chapter is, again, Reference [8]. The reader is encouraged 


to review this reference for further details. 


B. PROPELLER STRENGTH ANALYSIS--A HISTORICAL REVIEW 

Marine propeller blades present a special class of struc- 
tural problem. That problem lies in the difficulty of des- 
cribing a blade design in simple mathematical terms for 
subsequent analysis through various means. Until the "finite 


element method era", analytical methods, including the one 
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[Ref. 8] adapted for this study, relied heavily on practical 
experience of the propeller designer and semi-theoretical 
considerations. Analysis by these methods provided a cri- 
terion of stress rather than actual computation of stresses. 
These methods for predicting blade stresses were developed by 
using "beam" theory or "shell" theory. 

The use of elementary beam theory in propeller strength 
analysis was first adopted by Taylor [{Ref. 22]. He treated 
a blade as a cantilever beam attached to the propeller hub 
and loaded by thrust and torque forces distributed linearly 
over the propeller radius. His approach is often deemed a 
"modified beam theory" because he chose to calculate the 
direct stresses using the moment of inertia properties of 
expanded cylindrical blade sections with neutral axis parallel 
to the nose-tail (pitch-reference) line or chord of that 
expanded section. Reasonable estimates of stresses along the 
blade surface were achieved for the unraked, unskewed and 
narrow-bladed propellers of his time. 

As propellers "modernized" and became skewed and wider 
with increasing rake (mostly aft), modifications, improvements 
and alternatives to Taylor's theory were developed. Prin- 
cipally, modifications by Rosingh [Ref. 23] and Hancock 
{Ref. 24] proposed using moment of inertia properties of a 
blade section that was normal to the generating line of the 
axially projected blade outline. Romson [Ref. 25] later 


improved Taylor's theory for application to wide-bladed 
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propellers. Morgan [Ref. 26] provided an improved method 

for calculating the geometric properties of "modern" airfoil- 
shaped blade sections. Aernoldus and Keyser's [Ref. 27] 
“quasi-static” modeling of the propeller blade allowed 

for additional consideration to stresses induced by centri- 
fugal loading of the raked and skewed blade. The beam 

theory approaches to propeller blade stress analysis culminated, 
for all practical purposes, with Schoenherr's work [Ref. 8] 

in 1963. 

Alternatives to Taylor's beam theory approach, prior to 
1963, consisted of the application of "shell" theory to the 
propeller blade strength problem. This approach was first 
proposed by Conn [Ref. 28] and subsequently formulated by 
Cohen [Ref. 29] who modeled the blade as a helicoidal shell 
with variable thickness and infinite width. "Shell Theory” 
was utilized again in experimental studies by Connolly 


(Ref. 30] who, like his predecessors, was also forced into 


making an assumption about the behavior of the displacements 
of the blade sections (i.e., constant displacements normal 
to the constant pitch blade at each fractional radius dis- 
tance from the hub) beyond usual assumptions of shell theory. 
Essentially, his experimental results on one specific pro- 
peller contradicted the computational values. Attempts at 

a generalized numerical solution to Connolly's equations 
appeared in 1963 [Ref. 31] and 1964 (Ref. 32]. In 1968, 


Atkinson [Ref. 33], compared Connolly's results with currently 
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adopted cantilever beam methods and found, based on the 
inconsistency of results, that it was not possible to. 
recommend one method over the other in the blade strength 
design procedure; another approach was needed. 

Commenting on Atkinson's paper at that time, Sontvedt 
{Ref. 34] pointed out that, in view of these inconsistencies, 
the only approach to blade strength analysis which did not 
require very broad assumptions was the finite element tech~ 
nique. Developments in the method at that time were provid- 
ing a new and powerful tool for structural analysis. The 
propeller blade was just another application. Genalis [Ref. 
35] developed codes for the determination of displacements 
and stresses in a blade under hydrodynamic loads using the 
FEM technique and modeling the blade as a shell, a 3-D ele- 
ment mesh of tetrahedrons and rectangular prisms and, finally, 
a composite of shell and 3-D elements. As an aside, a 
finite "difference" solution to Connolly's analytical equa- 
tions was proposed in 1972 (Ref. 36]. In 1973, Atkinson 
{Ref. 37] reported the application of both hydrodynamic and 
centrifugal loads to a blade modeled by a thin-shell triangu- 
lar mesh and a thick-shell parabolic and cubic curved element 
mesh. The results of the triangular element were considered 
unsatisfactory. Another use of the thin-shell triangular 
element was reported by Sontvedt [Ref. 34] in 1974 using the 
SESAM-69 code (Ref. 38]. 

The need to model the blade correctly near the hub, 


where root stresses are usually critical, necessitated the 
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consideration of 3-D elements in lieu of thick/thin-shell 
elements (Refs. 39,40]. The use of the 4- and 10-noded 
tetrahedral element (i.e., TET-4, TET-10 respectively) to 
construct a blade mesh was conducted by Beek [Refs. 41,42]. 
He observed that improved accuracy of stress values, achieved 
by the use of these meshes, were overshadowed somewhat by 

the extensive storage capacity required for each analysis. 
Another "natural" improvement, from a geometrical standpoint, 
appeared in 1978. A general 3-D curved isoparametric element 
was incorporated in a computer code [Ref. 43] developed by 
Ma based on his previous formulation work [Ref. 44]. 

The finite element approach will continue to grow in use 
in propeller blade strength analysis with each successive 
improvement made to the basic elements which are used in the 
mesh generation of the blade. But, until the computer storage 
problem is resolved to the point where one mesh generation 
and subsequent stress analysis of one particular blade be- 
comes a minor processing task, the basic analytical techniques 
will continue to be a meritable "check" [Ref. 45] in the 
preliminary (or conceptual) phase of propeller selection/ 
design. In this context, the method formulated by Schoenherr 
and his colleagues twenty years ago is considered for adoption 
in determining a minimum required blade section maximum 


thickness-to-chord ratio. 
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C. SCHOENHERR'S METHOD 
1. Background 

Schoenherr's method is applicable to sreliwinaty 
(or, conceptual) propeller selection problems because it 
employs an assumed thrust and torque force loading distribu- 
tion for the propeller blade. This assumption is made at 
this stage of the ship's design because the exact wake 
velocity distribution at the propeller race is generally 
not known. 

Also, his method is applicable to propellers repre-~ 
sented by the Wageningen B-Screw Series because the blades of 
these propellers meet Schoenherr's criteria for the blade 
types covered by his formulation. Specifically, B-Screw 
Series blades have: 

1) a small constant rake angle of 15° over the entire 
blade radius; 

2) a constant pitch distribution over the propeller 
radius with the exception of the Z = 4 propeller whose pitch 
is slightly reduced near the hub; 

3} mild skew; 

4) linear distribution of blade section maximum thickness 
over the radius of the blade; 

5) aerfoil profile qualities where the nose-tail line 
or the chord of the blade section is approximately parallel 


to the pitch reference line. 
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2. The Blade Model 
Schoenherr models the propeller blade as a cantilever 

beam with unsymmetrical and variable area cross sections 
subjected to loading distributions of hydrodynamic and cen- 
trifugal forces. The following additional assumptions apply: 

1) Flexure theory applies. This subsequently implies the 
following: a) plane cross sections remain plane under load, 
b) Hooke's Law is valid, c) the blade material is homogeneous 
and isotropic, d) fibers are free to extend and contract 
independently of adjacent fibers, and e) stresses at a point 
arising from various forces superimpose. 

2) Shearing stresses and their effects are neglected. 
Only the direct stresses on a strength section are taken 
into account. 

3) The strength sections are taken to be the expanded 

1 cylindrical blade sections at various radial locations. 

4) The neutral axes of a strength section are straight 
lines passing through the centroid of the expanded cylindrical 
blade section and are parallel and normal, respectively, to 
the pitch reference line, and therefore, the chord, at each 
blade section. 

5) Bending Moments are applied in two planes which are 
mutually perpendicular to each other. One plane is normal 
to the pitch-reference line (and chord) of the strength 
section; the other is parallel. 

6) The angle between the principal axes of inertia and 


the neutral axes is zero. 
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Using this model and assumptions, Schoenherr applies 


the following formula for the evaluation of the direct fiber 


stress (To ,)) in a blade section at a radius r = ro? 


(M]} u {MJ w {F] 
gle. ae ee 8g ko 65 mey ) (6.1) 
. Tho To A 5 


(M} 0, = resultant bending moments in planes 
normal ((M]__) and parallel ((M] to 
the strength section's chord at ¥°= r 
(£t--Lbf) ; . 


[F_] = centrifugal force acting normal to the 
plane of the strength section at r=r 
and resulting from the centrifugal 
acceleration of the remaining blade 
element mass above that strength section 
(1lbf); 


Ww = coordinates of a point on the strength 
section's periphery with respect to 
that section's neutral axes system (L-n 
system) (ft); 


A(x_) = strength section's cross-sectional area 
m (£29) 


I = moment of inertia of the strength 
section with respect to the "2" axis 
(£t4) ; 


I = moment of inertia of the streng section 
with respect to the "n” axis (ft*); 


x = non-dimensional radius given by 


ae r/R- 


Since equation (6.1) indicates that the direct fiber stress 
is greatest at points on the periphery of the strength sec- 


tion, Schoenherr selects to examine four "critical points" 


81 


on the periphery where the fiber stress is likely to ce a 
maximum. These points are designated as (a), (2), G) and 
G) on Figure (6.1) and are specified by coordinates fate wig 
(u2,w2), (u3,w3) and (u4,w4) respectively in the "2-n" 
reference system. 


The values for [M] and IM] 9, are determined by 


no 
the following relation 


[M] 1 [Mplng + Men] no (6.2) 


[M],, = (Mp] (6.3) 


P° Lo + (M 


cb! 20 
where: 


= total bending moment due to hydrodynamic 
loading acting in a plane normal toa 
strength section's chord at r=r 
(£ftlbf) ; : 


[Mplo, = total bending moment due to hydrodynamic 
loading acting in a plane parallel to 
a strength section's chord at r=r 
(£tlb£); : 


total bending moment due to centrifugal 
loading acting in a plane normal to a 
strength section's chord at r=r 
(ftlb£) ; : 
total bending moment due to centrifugal 
loading acting in a plane parallel to 

a strength section's chord at r=r 
(ftlbf). 2 


3. Bending Moments Due to Hydrodynamic Loading 
The derivations of [MoJo and [My loy follow directly 


from Part I of Schoenherr's paper [Ref. 8: p. 83-89] and, 
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therefore, only key equations will be restated. References 
to equations which contain no decimal point apply to equations 
as numbered in his paper. 

Thrust and torque force are components of the hydro- 
dynamic "lifting" force acting on a blade. Using an assumed 
non-linear distribution of thrust along the blade radius 
given by equation (2), Schoenherr derives the following ex- 
pression for the bending moment due to thrust ((M,1,) which 


acts at a blade section located at ratius r = r_: 


° 
o5(x_) 
TR 2°00 
{M, ] = =: (6.4) 
t’o Z 9, (&) 
where: 
T = propeller thrust (lbf); 
R = propeller radius (ft); 
Z = no. of blades; 


functions of non-dimensional 
radius x evaluated at x_ = r/R 
and x, = 0.2 and given By eqia- 
tions (4) and (9). 


2 (%,) 04 (x) 


For the bending moment due to torque (IMJ) which 
acts at a blade section located at radius r = to? Schoenherr 


derives the following: 


IMJ, war a 1%) (6.5) 


Q, = propeller torque (ft-lbf); 


Z = no. of blades; 


V(X) 604 (%) 


functions of non-dimensional radius x 
evaluated at x_ = r_ /R and x, = 0.2 
and given by eQuati8ns (19) Yna (4). 


Figure (6.2) depicts the component resolution for 
[Mp] no and [Mp] 9, which results when equations (6.4) and 
(6.5) are imposed at a strength section at r = ry which has 
a pitch angle Bo: The following relations are derived as 


equations (42) and (43) in Schoenherr's paper: 


[Mp] no = [Mo], cos By + (Mol, sin By (6.6) 


[M1], sin By - [Mol cos By (6.7) 


4. Force and Bending Moments Due to Centrifugal Loading 


The derivation and expressions contained in this 
section constitute the author's modifications to the formu- 
lation in Part II of Schoenherr's paper. In Part II, Schoen- 


herr's derivations for [M and [M._,] are formulated for 


cb 
computation using a propeller drawing. This follows from 


eb) no Lo 


the fact that his method, which was funded by the American 
Bureau of Shipping, was intended to be used as that classi~ 
fication society's "designer's check" on adherence to the 

Bureau's strength criteria from a propeller blueprint. To 


evaluate the direct fiber stresses from equation (6.1) for 


the Wageningen B-Screw Series in accordance with Schoenherr's 
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method, (Mon) no and (Mop! g6 must be evaluated from expres- 
sions derived from the available information on blade sec- 
tion profiles and other geometric characteristics which 
are contained in Reference [2] and previously used in Chapter 
Vv. 

Consider Figure (6.3) where the centrifugal force 


{c of a blade element above a blade section at x = x_ = r/R 


Flo ° 
acts in a radial direction from the shaft centerline. Its 


line of action passes through point "N", which is on the same 


cylindrical surface as the blade section at x = x, = r/R, 


and through point "G", which is the center of gravity of the 
blade element above the blade section at x = x= r/R. 


From the figure, the follcwing expression is derived: 


(c (Cay cos G(x.) + (Chl, sin 6 (x,) (6.8) 


Fl 


[Ch], is shifted to point "N" and is decomposed into 
components [C,] cos ¢(x,) and [C,], sin c(x,)- The entire 
cylindrical surface in which the blade section at x = xs r/R 
and point "N" lie is now expanded into a flat plane for further 
consideration (see Figure (6.4)). In this configuration, 

[C1], cos (X,) is normal to this flat plane while [Ch], sin 6(x,) 
lies in this plane. 

Let point "O" be the location of the blade section's 
neutral axes system (i.e., the %£-n system). Then, the forces 
and moments due to the centrifugal reaction of the blade ele- 


ment above this section act at point "0" and are given by: 


87 


| : uNg 2UTOd uoTZOeSzaAUT pue ,D, APTAPAD Jo Ao_ZUA|D ¢°9 eanbrg 


88 


ANIT YOLVAaNAD ANIT FONIMGIA HOLId 


a/°rx = °x = x 3e uOT}OeS BpeTG e UO SuTpeoy TeurezurT f°9 ambtA 


89 


SNIT DJONTYaAIaa 


ANIT YOLWYAND 


[Fly = [Ch], cos ¢ (x,) (6.9) 
(Dol, = [C,], sins (x,) (6.10) 
[Mon], [Fod, + NO (6.11) 
(Mowlo = (Dol, * LO (6.12) 
where: 
[Fal = direct force, due to centrifugal action, 
acting on the blade section located at 
x = X, = r/R: 
(Dol, = shear force, due to centrifugal action, 
acting on the blade section located 
at x = es r/R: 
(Monle = bending moment at point "O" imposed by 
[Palo acting through point "N"; 
Mowso = torsional moment at point "0" imposed by 


[Duly acting through point "N". 


Since Schoenherr's method does not consider shear forces and 
their effects, [Do], and [Mowlo will not be considered in 
this modification. However, [Fo1, and (Mont, must now be 
computed for each blade section along the propeller's radius 
in order to account for their contributions to equations 
(6.2), (6.3) and, finally, in equation (6.1). 

The computation is derived as follows. Consider 
Figure (6.4). Again, (Fol, acts through point "N" and is 
normal (outward) to the plane of the figure. (Mo) is now 


° 
j resolved into components of the 2-n axes system as follows: 
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(Mop! no [Ch], cos C(x.) {py sinB,+q, cos B+ ¥.,) (6.13) 
[Man] eo = [C1], cos 5(x,) ta, sin 8 -p, cos Bo + Xoo} (6.14) 
where: 
Bo = pitch angle of the blade section at 
x = x = xr_/R; 
° ° 
x = distance of the blade sections centroid from 
| the generator line (ft); 
Yoo = distance of the blade section's centroid from 
the pitch reference line (ft); 
a distance to point "N" from point "P" 
parallel to the shaft axis at x = xX, = r/R 
(ft); 
a distance to point "N" from point "Pp" 


perpendicular to the shaft axis at 
x =X, = r/R (ft). 


The quantity 8 is found by the relation: 


i 
T 


tan Bs = . (P/D) , (6.15) 


| For the Wageningen B~Screw Series, (P/D) , is a constant along R 


except for propellers with 2 = 4. 


The quantity [Ch], is computed from the relation: 


7 wd V 


: Oo. 2. ee 3 
[C,], = 1728 SoG (2mn,) R* (%&Q), (6.16) 


91 


Ve = volume of the blade element above the blade 
section at x = x, = r/R (ft~); 

acg = 32.174 £t/sec*; 

wd = material weight density (1b£/in*) ; 

n, = propeller revolution rate (rps); 

(%)g = non-dimensional radial position of "G" 

for the blade element above the shaft 
axis; 

R = propeller radius (ft). 


At this point, only five quantities remain to be determined 
for the evaluation of the expressions of equations (6.10), 
(6.13) and (6.14). They are: 


1) cos ¢ (x,) 


2) Py ‘ 
' 3) Qo 
4) (X_), 
| 5) V5 


These quantities are determined by integration over the 


; blade element above the blade section, located at x = Xx, = r/R, 


from x = XG to x = 1.0. 
The values for Py and qo will vary with the location 
of "G" (from Figure (6.3)) which depends on x. Consider a 
radially thin slice of the blade element above the blade 
section, located at x = x, = r/R (see Figure (6.5)). This 
thin "slice" is located at a non-dimensional distance x from 
the shaft centerline where x, <x < 1.0. Figure (6.6) depicts 


this section expanded onto a plane. Let "g" be the centroid 


of that "slice". If Xg is the distance from the generator 
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line to "g" in feet and Yg is the distance from the pitch 
reference line in feet, then, from Figures (6.5) and (6.6), 


the following relations apply: 


i, = x Rtann (6.17) 

i,-a, = ig-%y sin By - Ve cos By (6.18) 
t = - si 19 

7 Xy cos 8, Yg sin By (6 ) 


where: 


n = rake angle at x; 


x _ = distance of "g" from point "P" parallel to 


g the shaft axis (ft); 
y = distance of "g" from point "P" perpendicular 
g to the shaft axis (ft). 


For the Wageningen B-Screw Series, rake angle n is a constant 
15° everywhere along the radius R. 

Now, to compute the volume of the blade element above 
the blade section located at x = x, = r/R, the following 
expression is used: 


v. = f RA(x) ax (6.20) 


To compute the non-dimensional radial position of 
"G" for the blade element above a blade section located at 


eS r/R, the following expression is used: 
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=X . 
Xo = > (6.21) 
if A(x) dx 
=x 
°o 


To compute the tangential position of "G" for the 


blade element above a blade section located at x = x, = r/R, 
the following expression is used: 
1 
f A(x) t_ ax 
x= g 
(6.22) 


And, finally, to compute the axial position of "G" 


for the blade element above a blade section located at x = Xe 
= r/R, the following expression is used: 
1 
A i -a 
L. (x) (Ag g) & 
oe iO 
Ay I (6.23) 
J A(x) dx 
*o 
Using the values just determined for (XG) 67 T, and 


A_, the following expressions are used to evaluate Py and 


a 


xX 
= (6.24) 


oa 


7 


qo = ASS x tan n (6.25) 


The expression for cos 5 (xX,) follows: 


P 


cos 5(x,) = *R (6.26) 


The formulation is now complete. Equations (6.10), 
(6.13) and (6.14) can now be evaluated for any blade section 
located at x = x, = r/R. From here, equations (6.2) and 
(6.3) are evaluated. Finally, using the results from these 
equations and equation (6.9), equation (6.1) can be evaluated 
for the four points, specified by coordinates (ul,wl), 
(u2,w2), (u3,w3) and (u4,w4), at any location x = x, = ro/R- 


From the development discussed in Chapter V, the 


values for A(x), x_, y_, I 


g g ae and Ip, are readily available. 


D. ALGORITHM FOR THE CONSTRAINT 
1. Theory 
Schoenherr's formulation with the modifications just 

derived can be used in determining the minimum required 
equivalent blade section maximum thickness~-to-chord ratio 

* s ry r+ . 
((t*/c) a5p min) for use in the constraint G,,(X) < 0 given 
by equation (4.8). The procedure employed is as follows: 

pelea Fe , 
1) assume an initial value for (t /C) a5p min U8ing 


equation (3.21); 


2) increase this value by a small amount; 


97 


3) using (t*/c) s6p ae obtained from step (2), generate 
a distribution of minimum required blade section thicknesses 
(tein? for blade sections at specified points along the 
propeller radius, say at r = .2R, .3R, .4R, .5R, .6R, .7R, 
-8R and .9R; 

4) determine all blade section properties to include: 
a) cross-sectional area, b) chord length, c) centroid location, 
d) moments of inertia with respect to the principal axes 
system (i.e., &£-n system), e) coordinate values for the four 
critical points defined in the previous section; 

5) compute the hydrodynamic bending moment components 
[Mp] no and [Mplog at radius locations just specified; 

6) compute the values of the centrifugal force [Foly and 


the bending moment components Mob!) no and Mop! acting on 


Lo 
blade sections at radius locations just specified; 
7) calculate the direct fiber stresses at all four 
critical points for all radius locations specified in step 
(3); 
8) check the following condition on the calculated fiber 


stress at all four critical points at all specified radius 


locations using: 


lo], < 144° 8 (6.27) 


9) if the maximum allowable stress (S.) for the material 
is exceeded, then return to step (2) and repeat. Otherwise, 


proceed to next step. 
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10) since the minimum required equivalent blade section 
Maximum thickness-to-chord ratio assumed in step (2) has 
produced blade sections of adequate strength, evaluate the 
constraint given by equation (4.8). 

2. Coding Details 
The algorithm just outlined is incorporated into the 
main FORTRAN subprogram SUBROUTINE STRCNK. This subprogram, 
in turn, executes the algorithm through sequential calls to 
other key FORTRAN subprograms. These subprograms are listed 
as follows: 

1) SUBROUTINE TDIST--accomplishes step (3); generates, 
at specified radius values, a distribution of minimum required 
blade section maximum thicknesses (trim for the assumed 
value of (t*/c) s6p min? 

2) SUBROUTINE BLDPRP--accomplishes step (4); described 
previously in Chapter V; 

3) SUBROUTINE HYDLD--accomplishes step (5); computes the 
hydrodynamic bending moment components, given by equations 
(6.6) and (6.7), at specified radius locations; 

4) SUBROUTINE CNFGLD--accomplishes step (6); computes the 
centrifugal force and bending moments, given by equations 
(6.9), (6.13) and (6.14) respectively, at specified radius 
locations; 

5) SUBROUTINE SIGNDS~-accomplishes step (7); computes 
direct fiber stresses, given by equation (6.1), for all four 


critical points at every specified radius location. 
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During the remaining steps of SUBROUTINE STRCNK, the 
condition on allowable stress, given by equation (6.27), 
is checked at all critical points of blade sections located 
at specified radius locations (again, r = .2R, .3R, .4R, .5R, 
-6R, .7R, .8R and .9R). The final calculation made is 
that for the constraint given by equation (4.8). 

Again, extensive use of common blocks, for passing 


data from one subprogram to another, is apparent upon examina- 


tion of the codes just cited. Comment cards are used 
throughout. 
E. SUMMARY 

b The end of this chapter marks the completion of all 


prerequisite background and formulation discussions on the 
application of COPES/CONMIN to propeller selection problems 
involving the Wageningen B-Screw Series. From this point, 


each specific Design Case can now be solved as an optimiza- 


tion problem. 
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VII. DESIGN CASE NO. 1--PROGRAMMING AND COMPARISONS 


A. INTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
propeller selection problems which use the "thrust" approach. 
First, the thrust approach to the propeller selection problem 
is formulated. Then, a review of a previous author's solution 
to this problem is presented. Four variations to this pro- 
peller selection problem are solved by COPES/CONMIN. The 
chapter is completed with a presentation and discussion of 


the results from the four variations. 


B. THRUST APPROACH FORMULATION 
1. Design vector XI 
As previously pointed out at the conclusion of Chap- 
ter III, Design Case No. 1 constitutes a propeller selection 


problem which is solved by the thrust approach. In this ap- 


proach, the effective horsepower (Pp) and the ship's speed 

(V) are specified by the designer. From the viewpoint of 
optimization, the quantities Pr and V become preassigned 
parameters. This reduces the design vector X (see Figure 4.1) 


tos 


(7.1) 


I 
) 
» 
fos] 
> 
ie) 
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Having specified Pe and V, all of the design variables, 
as listed in equation (7.1), are not independent. Recalling 


equations (3.3) and (3.13), the following relationship 


results: 
oe ae CIR we ee aes 
E (i-wt) "R° "9 ° B56 ° 80 ; 


Rearranging terms, this equation becomes: 


P 


- {izwt) . "ez . 350-60 


Considering that the open water efficiency (ng) is evaluated 
prior to the computation of (=n) + or OBJ; 2° then both Np 
and Q,g are not independent design variables. One must 

be selected as the independent design variable. Then, 

the other variable becomes dependent on the one just 
selected. 


For this study, N, is selected as the independent 


P 
design variable. This choice will reduce the design vector 
X for propeller selection problems using the thrust approach 


to the following: 


(7.4) 


>] 
) 
re 
> 
1°) 
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Finally, equation (3.8) implies an alternative defini- 
tion of X as given in equation (7.4). The design vector for 
Design Case No. 1 propeller selection problems is, therefore, 


defined as: 


xl = A,/Ay (7.5) 


2. Powering Constraint 


Having determined the design vector Xl, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 

restriction constitutes the remaining constraint Gy 5 (X) men- 


tioned in Chapter IV. 
Simply stated, the selected propeller, as defined by 


XI, must develop enough thrust (T) so that the powering require- 
ment, specified by Pp and V, is met. Using equation (3.2), 
the thrust developed by the propeller can be specified in terms 


of thrust horsepower (Pin) as: 


_ Tf V(l-wt) 
(Po) gey ~ “S550 (7.6) 


From equation (3.9), it follows that: 


2 4 
pn, D Ky 


(Pn) gey aera <7, amelie V(l-wt) . (7.7) 


Using equation (3.12), the developed thrust horsepower can 
be defined in terms of developed effective horsepower given 


by: 


te . (Pa) aey (7.8) 


(PE) aev 
The restriction imposed by the thrust approach method, 

where Py and V are specified, can now be stated as: 

(7.9) 


Ph oS (Px) 


E dev 


Rearranging equation (7.9), the constraint Gy, (x1) follows: 


(P.) 
- peer < 9 (7.10) 


G(Xl)°° = 1 
12 = 


With the design vector XI and G, 9 (XT) defined, the 
propeller selection problem represented by Design Case No. l 


can be stated under one equation as: 


Max; mize: F(xXI) = OBJ 
Subject to: G, (x1) < 0 5 eae Le (7x1) 
lower upper . _ 
x1; < x1, < x1; 1. Lys snes 
C. PREVIOUS SOLUTIONS 
Triantafyllou [Refs. 3,21] considered a propeller selection 
problem represented by Design Case No. 1. In his example 


problem, the following parameters were specified: 
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1) v = 1.139 x 107© (m?/sec) = 1.22613 x 107° (£t?/sec) 
2) wt = .22 

3) ta = .19 

= 1.025 


5) noscrw = l 


6) Z=5 

7) Py, = 18153 (hp) 
8) v= 24 (knots) 
9) D= 22.0 (ft) 


10) Az/Ay = .85 
The hull under study in his example had the following dimen- 
sions: 
1) Length = 710 (£t) 
2) Draft = 30 (ft) 
3) Beam = 100 (ft) 
For his analysis, the design vector contained two varia- 


bles and was specified as: 


Using an iterative scheme [Ref. 21: p. 79] to solve two equa- 
tions in two unknowns, he maximized the open water efficiency 


(n,) to obtain the following results: 


P/D = 1.1651 
N, = 104 (rpm) 
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P = 25544 (hp) 


ae -6676 


For future comparisons, equations (3.8) and (3.3) give: 


J 8286 


1290000.0 (ft-lbf) 


Qs 
Triantafyllou's results are summarized in Table (IV). 


D. SOLUTIONS BY COPES/CONMIN 
The propeller selection problem, as stated by equation 
(7.11), is now solved by COPES/CONMIN. Four solution varia- 
tions are considered. 
The first and second variations attempt to reproduce the 
solution given by Triantafyllou. The design vector XT 
(NDV = 2) is used in both cases. One variation uses SUBROUTINE 
STRCNA to evaluate the constraint Gy (XT) given by equation 
(4.8). The other uses SUBROUTINE STRCNK to determine G,. (XT). 
The remaining two variations will solve the propeller 
selection problem using the design vector XI (NDV = 5) defined 
in equation (7.5). Again, one variation uses SUBROUTINE 
STRCNA; the other, SUBROUTINE STRCNK. 
In all variations, the following parameters are used: 
1) Temp = 59 (°F) 
2) p = 1.9384 (1bf-sec*/£t*) 
3) v = 1.2285 10 (£t*/sec) 
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4) .247 (psia) 


Pwatvap ~ 


5) = 14.7 (psia} 


Patm 
6) wt = .22 


7) td = .19 


8) = 1.025 


NR 
9) noscrw = l 


10) 4h = 19.0 (ft) 
cy 

11) Diim = 22.0 (ft) 

12) 225 


13) promat = 5 (stainless steel; see Table (II)) 
14) P, = 18153 (hp) 
15) V = 24.0 (knots) 
All of the above are initialized in the input phase (ICALC = l) 
of each SUBROUTINE ANALIZ pertaining to each variation. 


The constraint Gy 4 (XT) or Gy > (XT) is evaluated by SUBROU-~ 


TINE BLPOW1 which appears in the execution section of each 
SUBROUTINE ANALIZ. 
1. Variation 1 
a. Programming Details 
Since this variation uses the design vector XT, 
the following design variables of XI become paransters and 
are specified in the input section of SUBROUTINE ANALIZ (ICALC = 1) 
as: 
1) D = 22.0 (£t) 
2) A,/AG = .85 
3) (t*/c) a5p = .0348 (from equation 3.21). 
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Re 


For constraints, the following are used: 
G; (XT) < 0 j =1,...,8, 12 


Only nine of twelve constraints are evaluated (NCON = 9). 
Obviously, some of the twelve constraints are redundant since 


D, Ap/Ags and (t*/c) 45p have been specified. 


lower 


PEpOr; and lower (XTY 


The upper (XT; ) limits on 


the design variables J and P/D are set to be: 


1.6 


° 
a 
[A 
Cy 
JA 


-4 < P/D < 1.4 


These upper and lower limits are specified in the COPES 
control card deck on card image F under respective fields 
VUB and VLB. The initial value for each design variable 
(XT; ) is also assigned on card image F under the field labeled 
X. The first list of card images in Appendix D lists all of 
the COPES control cards used for this variation and variation 
2. These cards also specify the locations of the design 
variables in the common block GLOBCM (see Table (III)) as 
well as the locations of the constraints and their boundaries. 
Further details on the COPES control card requirements and 
the format of each card are contained in Reference [7]. 

An examination of SUBROUTINE ANALIZ for this varia- 
tion, found in Appendix C, shows the calling statement made 


to SUBROUTINE STRCNA. 
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b. Results 
The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed first in Appendix E. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (IV). 
2. Variation 2 
a. Programming Details 
Everything discussed above for the first variation 
applies here with one exception. An examination of SUBROUTINE 
ANALIZ for the second variation, found in Appendix C, shows 
the calling statement made to SUBROUTINE STRCNK. 
b. Results 
The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed second in Appendix E. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (IV). 
3. Variation 3 
a. Programming Details 
This variation uses the design vector XI. For 


constraints, the following are used: 
G, (x1) C0. Soe ipwag eh? 


All twelve constraints are evaluated (NCON = 12). 


lower 


PPET) and lower (Xl; 


The upper (x1; ) limits 
on the design variables D, P/D, Ap/Ao: (t*/c) 75R! and J 


are set as: 
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These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (X15) 

is also assigned on card image F under the field labeled xX. 
The second list of card images in Appendix D lists all of 

the COPES control cards used for this variation and variation 
4. These cards also specify the locations of the design 
variables in the common block GLOBCM (see Table (III)) as 
well as the locations of the constraints and their boundaries. 
Further details on the COPES control card requirements and 
the format of each card are contained in Reference [7]. 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix C, shows the calling statement 
made to SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed third in Appendix E. Re- 
sults for this variation of thepropeller selection problem 


are tabulated in Table (IV). 
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4. Variation 4 
a. Programming Details 
Everything discussed above for the third varia- 
tion applies here with one exception. An examination of 
SUBROUTINE ANALIZ for the fourth variation, found in Appendix 
C, shows the calling statement made to SUBROUTINE STRCNK 
instead of SUBROUTINE STRCNA, 
b. Results 
The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed last in Appendix E. Re- 
sults for this variation of the propeller selection problem 


are tabulated in Table (IV). 


E. DISCUSSION 

Overall, the results achieved in all variations compare 
reasonabiy well to the solution obtained by Triantafyllou. 
However, the following points can be made. 

Variations 1 and 2 give the same results. This was ex- 
pected in view of the fact that, even though constraints 
Gg (XT) through G, 4 (XT) were evaluated, these constraints were 
not considered in the optimization search conducted by CONMIN. 

In variations 3 and 4, the diameter (D) was driven to 
the limit (Di im): This bears out a fundamental rule in pro- 
peller design, i.e., the larger the propeller diameter (D), 
the greater the open water efficiency (n,)- 

The minimum required equivalent blade section maximum 


thickness-to-chord ratio ((t*/c) a5R min)? computed in 


lll 


variation 3, is substantially smaller than the one computed 
for variation 4. As pointed out in Reference (2], the empiri- 
cal relation, expressed by equation (4.9) and derived from 
equation (70.x) (Ref. 46: p. 620], does not take into account 
the effects of centrifugal loading. These effects include, 
specifically, the direct stresses imposed by the inertia load 
of the blade and the bending moments which result from rake 
and skew of the blade. Therefore, the algorithm developed in 
Chapter VI should, and does, produce a larger value for 
(t*/c) Jo5p- 
A final observation on the results concerns the values 
of the open water efficiency. The "optimum" open water effi- 
ciency (n,) achieved by Triantafyllou is lower than those 
achieved in variations 1 and 2. A possible reason for this 
might be the neglection of the term "dRe/dJ" in Triantafyllou's 
formulation of the analytical expressions [Ref. 21: p. 71] 
that he used in his analysis. The difference in the open 
water efficiencies subsequently accounts for the differences 
in the propeller revolution rate (N,) and the delivered 


torque (Q,) when the relation in equation (3.3) is considered. 
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TABLE IV 


Design Case No. 1--Results 


GROUP I'TEM TRIANTA~ 
P. 18153.0 18153.0 18153.0 18153.0 18153.0 
Given Vv 24.0 24.0 24.0 24.0 24.0 
Design D 22.0 22.0 22.0 
carga Ap/Bp .85 .85 85 
fied (EN/e). oe .0348 .0348 
D 21.9991 21.9659 
Design P/D 1.1651 1.0036 1.0036 .9981 1.0071 
Variables Ae/Ay 8205 —»-«. 8149 
(t*/c) Jen .0330 0642 
3 .8286 =. 7371 7371 .7343 .7394 
Maximize n,, .6676  .7091 7091 .7109 7066 
Di im (22.0) (22.0) 22.0 22.0 
Restric~ AYA, (0.5258) (0.5258) 5269 5267 
tions (t*/e) sent (0.21266) (0.50761) .021998  .053259 
N, 104 116.9 116.9 117.4 116.7 
Other Q, 1290000 1080451. 1080451 1064574 1084003 
Py 25544.0 24051.3  24051.3 24030.3 24094.8 
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VIII. DESIGN CASE NO. 2--~PROGRAMMING AND COMPARISONS 


A. INTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
propeller selection problems which use the "power" approach. 
First, the power approach to the propeller selection problem 
is formulated. Then, a review of a previous author's solu~ 
tion to this problem is presented. Four variations to this 
propeller selection problem are solved by COPES/CONMIN. The 
chapter is completed with a presentation and discussion of 


the results from the four variations. 


B. POWER APPROACH FORMULATION 
1. Design Vector XZ 
As previously pointed out at the conclusion of 

Chapter III, Design Case No. 2 constitutes a propeller selec- 
tion problem which is solved by the power approach. In this 
approach, the delivered torque (Q,) and the propeller revolu- 
tion rate (N,) are specified by the designer. From the view- 
point of optimization, the quantities Qs and Np become pre- 
assigned parameters. This reduces the design vector X (see 


Figure (4.1)) to: 


*| 
N 
o 


(8.1) 


Having specified Q, and Nye all of the design varia- 
bles, as listed in equation (8.1), are not independent. Re- 
calling equations (3.3), (3.11) and (3.13), the following 


relationship results: 


P, = (i-td) - ot - anos . Np (8.2) 
E (i-wt) ° "R 21K 550 ° 60 : 


Rearranging terms, this equation becomes: 


Q N 
= ——aoe 8 © acne 26 ee S ° P (8 3) 
Vv (l-wt) R on, D K 550 60 . 


Considering the relations for Ke and Ky in equation (3.17), = 
then both PE and V are not independent design variables. 
One must be selected as independent, while the other becomes 
dependent on the one selected. 

For this study, V is selected as the independent de- 
sign variable, This choice will reduce the design vector xX 
for propeller selection problems using the power approach to 


the following: 


xX = P/D (8.4) 
Ap/A 
* 
(ENS) 955 
Finally, equation (3.8) implies an alternative defini- 


tion of X as given in equation (8.4). The design vector for 
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Design Case No. 2 propeller selection problems is, therefore, 


defined as: 


V 
z 

P/D (8.5) 
Ag/Ay 

a Aa ie | 


iI 


2% Powering Constraint 


Having determined the design vector X2, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 
restriction constitutes the remaining constraint G(X) men- 
tioned in Chapter IV. 

Simply stated, the selected propeller, as defined by 
X2, must absorb at least all of the power delivered to it 


(Pp) which is specified in terms of Qs and N Using equation 


P* 
(3.3), the power absorbed by the propeller can be specified 


in terms of delivered horsepower (Pp) as: 


(P 


N. 
= Pam 
p) absorb - 550 60 (8.6) 


From equation (3.10), it follows that: 


2 5 
Ky fe) Rp D 2m Np 


(Pp) absorb - 550 ° 60 (8.7) 


But, equation (3.3) also defines the power delivered to the 


propeller as: 
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P 
D 1 ns) (8.8) 


The restriction imposed by the power approach method, 


where Qs and N, are specified, can now be stated as: 


P 


Py & 'P 5) abscess (8.9) 


Rearranging equation (8.9), the constraint Gy 9 (X2) follows: 


_ (P,) 
(XZ) = 1 ~ —Biabsorb < 9 (8.10) 


G 
12 D 


Further simplification of equation (8.10) gives: 


Q 
Gothen -S, doses ei ng (8.11) 
12 Qo = 
With the design vector XZ and Go (X2) defined, the 


propeller selection problem represented by Design Case No. 2 


can be stated under one equation as: 


Minimize: F(X2) = OBJ 


ize 
Subject to: G, (X2) oD) fee lgdew ke (8.12) 
lower upper —_ 
X25 < x2, <2; i aS 


C. PREVIOUS SOLUTIONS 
Markussen [Ref. 4] considered a propeller selection prob- 


lem represented by Design Case No. 2. In his example problem, 


the following parameters were specified: 
1) Temp = 18 (°C) = 64.4 (°F) 


ade Pyatvap 


3) Paotm = 1-01312856 (bars) = 14.6974 (psia) 


= 0.0206411 (bars) = .29943921 (psia) 


4) noscrw = 1 


5) Hog = 6.7 (meters) = 21.9827 (£t) 


6) Z= 6 

7) Py = 18.9 (MegaWatts) = 25344.9 (hp) 
8) Na = 110 (rpm) 

9) v4 15.65 (knots) 


For his analysis, the design vector contained three 


variables and was specified as: 


| 


P/D 


AB/AG 


A restriction for the minimum required expanded area ratio 
((Ap/AQ) min) * given by equation (4.3), was also considered. 
This imposed a constraint given by equation (4.4). 

Using an iterative scheme [Ref. 4: p. 110] to solve three 
equations in three unknowns, Markussen maximized the open 


water efficiency i? to obtain the following results: 


J = .61095 
P/D = .864380 

2 je 
Ap/Ap = 36.1861/40.6123 (m°/m") 
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= .891012 


No = -654391 


For future comparisons, equations (3.8) and (3.3) give: 


D 


7.19091 (meters) = 23.593375 (ft) 


Q 1210130.0 (ft-1bf) 


Ss 


Markussen's results are summarized in Table (V). 


D. SOLUTIONS BY COPES/CONMIN 

The propeller selection problem, as stated by equation 
(8.12), is now solved by COPES/CONMIN. Four solution varia- 
tions are considered. 

The first and second variations attempt to reproduce the 
solution given by Markussen. The design vector XM (NDV = 3) 
is used in both cases. One variation uses SUBROUTINE STRCNA 
to evaluate the constraint G, 2 (XM) given by equation (4.8). 
The other uses SUBROUTINE STRCNK to determine Gy > (XM) . 

The remaining two variations will solve the propeller 
selection problem using the design vector X2 (NDV = 5) de- 
fined in equation (8.5). Again, one variation uses SUBROUTINE 
STRCNA; the other, SUBROUTINE STRCNK. 

In all variations, the following parameters are used: 

1) Temp = 64.4 (°F) 
2) p = 1.9892 (1bf-sec*/et‘) 


5 


3) v= 1.1900 x107> (£t*/sec) 


4) = .2994 (psia) 


Pwatvap 


FR eee etn Ng Ly 


— 7 7@ yn 0. 


5) = 14.697 (psia) 


Patm 
6) wt = .22 
7) ta = .19 


8) Vs ™ 1.025 


9) noscrw = 1 


10) Nog = 21.9827 (ft) 
11) Di im = 30.0 (ft) 
12) 2=6 


13) promat = 5 (stainless steel; see Table (II)) 


1210130 (ft-1lbf) 


14) Q, 


15) N 110 (rpm) 


P 
All of the above are initialized in the input phase 


(ICALC = 1) of each SUBROUTINE ANALIZ pertaining to each 
variation. 

The constraint G) 9 (X2) or G,,(XM) is evaluated by 
SUBROUTINE BLPOW2 which appears in the execution section of 
each SUBROUTINE ANALIZ. 

l. Variation 1 

a. Programming Details 
Since this variation uses the design vector XM, 
the following design variables of X2 become parameters and 
are specified in the input section of SUBROUTINE ANALIZ 


(ICALC = 1). The ship's speed (V) is specified as: 


Vv 


V,/ (i-wt) 


15.65/(1 =- .22) 


20.0641 (knots) 
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Markussen elected to use the standard Wageningen blade sec- 
tion maximum thickness-to~chord ratios. Since the egateavent 
t/c is given as a function of Z and Ap/Ap (see equation (3.21)), 
then (t*/c) sep is calculated during each analysis (ICALC = 2) 


by the following relation: 


(t*/c) gen = (t/C) 5p 


For constraints, the following are used: 


G, (XT) < 0 ji 2 Leng, 105. 12 


Only ten of twelve constraints are considered (NCON = 10). 
Constraints Gy (XM) and G1 (™) are redundant since no limit 


on the propeller diameter (D ) appears as a parameter in 


lim 
Markussen's formulation and (t*/c) 75R Was taken to be the 


Wageningen standard. 


lower 


The upper (xy PPS) and lower (XM; ) limits 


on the design variables J, P/D and Ap/Ay are set to be: 


Ol < J < l.l 
24 <  P/D < 1.4 
4 < ASAD < lel 


These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 


and VIB. The initial value for each design variable (XM; ) is 
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also assigned on card image F under the field labeled x. 

The first list of card images in Appendix G lists all of the 
COPES control cards used for this variation and variation 2. 
These cards also specify the locations of the design varia- 
bles in the common block GLOBCM (see Table (III)) as well as 
the locations of the constraints and their boundaries. Further 
details on the COPES control cacd requirements and the format 
of each card are contained in Reference [7]. 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix F, shows the calling statement 
made to SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, 
performed by COPES/CONMIN, is listed first in Appendix H. 
Results for this variation of the propeller selection problem 
are tabulated in Table (V). 

2. Variation 2 
a. Programming Details 

Everything discussed above for the first variation 
applies here with one exception. An examination of SUBROUTINE 
ANALIZ for the second variation, found in Appendix F, shows 
the calling statement made to SUBROUTINE STRCNK instead of 
SUBROUTINE STRCNA. 

b. Results 
The output from the optimization/analysis, per- 


formed by COPES/CONMIN, is listed second in Appendix H. 
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Results for this variation of the propeller selection problem 


are tabulated in Table (Vv). 
3. Variation 3 
a. Programming Details 
This variation uses the design vector X2. For 
constraints, the following are used: 


G, (X2) < 0 jim detec 2 


All twelve constraints are evaluated (NCON = 12). 


lower 
i 


upper 
i 


on the design variables V, P/D, Ap/Ay: (t*/c) 75R! and J are 


The upper (X2 ) and lower (X2 ) limits 


set as: 


10.0 < V < 100.0 (ft/sec) 
40< P/D < 1.4 
4 0< Ap/Ag < Lal 

.003 < (t*/c) 75R < 50 
-O1 < J < 1.1 


These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (X25) 

is also assigned on card image F under the field labeled X. 
The second list of card images in Appendix G lists all of the 
COPES control cards used for this variation and variation 4. 


These cards also specify the locations of the design variables 
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in the common block GLOBCM (see Table (III)) as well as the 
locations of the constraints and their boundaries. Further 
details on the COPES control card requirements and the format 
of each card are contained in Reference [7]. 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix F, shows the calling statement 
made to SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed third in Appendix H. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (V). 

4. Variation 4 
a. Programming Details 

Everything discussed above for the third varia- 
tion applies here with one exception. An examination of SUB- 
ROUTINE ANALIZ for the fourth variation, found in Appendix 
F, shows the calling statement made to SUBROUTINE STRCNK 
instead of SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed last in Appendix H. Results 
for this variation of the propeller selection problem are 


tabulated in Table (V). 


E. DISCUSSION 
The results achieved in variations 1 and 2 compare ex- 


tremely well to the solution obtained by Markussen. As 
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See 


pointed out in the discussion in Chapter VII, variations 1 

and 2 are expected to give the same results for the vector 

XM. Obviously, the values obtained for J and P/D, as well 

as those for D, and Rn* eR are very close to the values 
generated in Markussen's example. However, the values for 
Ap/Ay are somewhat different. It is interesting to note that 
the value obtained in variations 1 and 2 (and, for that matter, 
variations 3 and 4) is, essentially, the limiting value for 
Ap/Ay: as given in Table (I), for 2 = 6. Markussen's value 
for Ap/Ay (i.e., .891012) exceeds the limit (i.e., .80) in 
this table. 

As pointed out at the end of Chapter VII, the minimum 
required blade section maximum thickness-to-chord ratio 
((t*/¢) o6p min)? » computed in variations l and 3, is sub- 
stantially smaller than the one computed for variations 2 
and 4. Again, the same explanation applies here as well. 

The results of variations 3 and 4 differ somewhat from 
Markussen's results. The reason for this is simply that Va 
(or V) has not been specified as a parameter. Consequently, 
a higher value for the advance ratio (J), which corresponds 
to a higher open water efficiency (node has been found in 
the optimization search. This result can be interpreted 
in the following way. Given: 

1) a six-bladed Wageningen propeller (Z = 6) which is made 
out of stainless steel (promat = 5); 
2) a power train delivering 25344.9 (hp) at a rate of 


110 (rpm); 
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3) a hull with a wake fraction (wt) equal to .22, a 


thrust deduction (td) equal to .19 and a shaft centerline 
depth (h.,) of 21.98 (ft), . 

then, the selected propeller, as defined by X2, can drive 
this hull at a maximum speed of V when the hull has a maximum 


resistance given by Poe 
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Design 
Variables 


Maximize 


Restrictions Ap/Ay min 


Other 


TABLE V 


Design Case No. 2--Results 


ITEM MARKUSSEN 
Py 25344.9 
Q, 1210130 
N, 110 
Vp 15.65 
Vv 
J .61095 
V 
Va 
P/D . 864380 
Ae/Ay .891012 
(t*/¢) 75R 
1p .654391 
DL im 

(574729 
(t*/C) 5p ; 
D 23.53 
*E - 
Rn* 7. 6.47810 


VARIATIONS 
" 2 3 
25344.9 25344.9  25344.9 
1210130 1210130 1210130 
110 110 110 
15.65 15.65 
20.0641 20.0641 
6475 6475 9927 
30.9138 
24.1127 
.9036 9036 1.1986 
.8018 8018 .7946 
.0397 .0397 .0499 
6660 6660 7643 
(30.0) (30.0) 30.0 
.5070 .5070 4622 
(.02729) (0647)  .02706 
22.25 22.25 22.36 
(14057.3) (14057.3) 19945.5 
5x10’ 5x20’ 5x10’ 
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25344 .9 
1210130 
110 


.8753 
29.5355 
23.0377 

1.0308 

. 7986 

.0638 


. 7330 


30.0 
- 4266 
-0638 


24.23 
20653.7 
5x 107 


IX. DESIGN CASE NO. 3--PROGRAMMING AND COMPARISONS 


A. INTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
a propeller selection problem where "matching" is desired. 
First, the “matching” approach to the propeller selection 
problem is formulated. Then, a review of a previous author's 
solution is presented. One variation to this propeller 
selection problem is solved by COPES/CONMIN. The chapter is 


completed with a presentation and discussion of the results. 


B. "MATCHING" FORMULATION 

1. Design Vector x3 

Design Case No. 3, the final powering problem con- 

sidered in this study, constitutes a propeller selection prob- 
lem solved by the "matching" approach. In this approach, 
the hull's effective horsepower (Pp) and speed (V), the 
delivered torque (Q,) and the propeller revolution rate 
(Np) are specified by the designer. This reduces the design 


vector X (see Figure (4.1)) to: 


=p 
P/D 
A, /A 
* 
(t*/e) a5R 


red 
" 


(9.1) 


For this study, the design vector % is reduced further 


by eliminating the propeller diameter (D) as a design variable. 
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That is, D will also be specified by the designer so that 


the design vector for Design Case No. 3 is defined as: 


P/D 
3 o= (Ava, (9.2) 
(A/S) 78K 


2. Powering Constraint(s) 


Having determined the design vector X3, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 
restriction constitutes the remaining constraint Gy 9 (X) men- 
tioned in Chapter IV as well as an additional constraint. 

In the "matching" problem, the selected propeller, 
as defined by X3, must satisfy two conditions. First, it 
must develop, as a minimum, the effective horsepower (Pp) 
as imposed by the design specification. Citing the formula- 
tion previously derived in Chapter VII, this condition can 


be stated as: 


Pe < (PR) Gey (9.3) 
The constraint Gy 9 (X3) follows accordingly as: 
(P,) 
Gis ikay Si ive Bey eG (9.4) 
12 PE —- 


For the second condition, the selected propeller can 


only absorb, as a maximum, the delivered power (P)) as 
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specified by the designer. The formulation is the same as 
that in Chapter VIII except that the inequality signs are 
reversed. The condition is stated as: 

(P <P 


(9.5) 


D) absorb D 


In defining a constraint G1 3(X3), another location, say 
location 24, in the GLOBCM block (see Table (III)) would be 
assigned. But, considering the fact that constraint Gy (X) 
will not be used because the propeller diameter (D) is speci- 
fied, there is no reason why Gy (X3) cannot be redefined, for 


this Design Case only, as: 


(P.) 
_Diabsorb _ |; . g (9.6) 


P = 


G,(X3) = 
ss D 


Further simplification of equation (9.6) gives: 
ay = P 
Gg(X3) = 5-- 1 < 0 (9.7) 


In reality, the constraints just defined should be 
equality constraints. The word "match" does infer equality 
in some sense. However, as previously stated in Chapter II, 
the version of COPES/CONMIN used in this study does not 
directly handle equality constraints. But, since CONMIN 
attempts to minimize constraints in the optimization search, 


it will be assumed that a "match" can be achieved. 
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With the design vector X3 and the constraints 


G) > (X3) and Gy (X3) defined, the propeller selection problem 


represented by Design Case No. 3 can now be stated under 


one equation as: 


Minimize: F(X3) = OBJ, 
Subject to: G,(K3) < 0 j= 1,...,12 (9.8) 
as ek, Sxaunee Abe gs 


C. PREVIOUS SOLUTIONS 


The propeller selection problem considered by Vassilopoulos 


[Ref. 18] actually represents a propeller "design" problem 


using the "power" approach. The example problem which he 


elected to solve is taken from that posed by the International 


Towing Tank Conference (ITTC) Propeller Committee. This 


problem is concerned with the determination of propeller 


thrust (T), diameter (D) and speed of advance (Vy) (or ship 


speed (V)) for a single-screw cargo ship where: 


1) 


2) 
3) 
4) 
>) 


power available to the propeller (i.e., Py) is 
30,000 (hp) 


Z2= 6 
105--110 (rpm) 


N 
Diim = 23 (£t) 
Roe = 19 (£t) 


The variation of ship speed (V) and of hull effective power 


(PE) ‘ 


thrust deduction factor (td) and the wake fraction 


(wt) is also given (Ref. 18: p. 20]. 
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The results from Vassilopoulos' propeller design exercise 
produced a propeller that is "matched" at the following values: 
1) P,Q = 21292.6 (hp) 
2) V = 24.24 (knots) 


3) Q< = 1500606.75 (ft-1bf) 
4) Np = 105 (rpm) 
5) Py = 30000.0 (hp) 


His propeller "design" was based upon the following specified 
parameters: 

1) Temp = 59 (°F) 

2) p = 1.9905 (1bf-sec*/st*) 


3) .247 (psia) 


Pywatvap a 
4) wt = .22 
5) td = .1725 
6) noscrw = 1 
7) he = 19.0 (ft) 

8) Z2=6 

9) promat = 5 (stainless teel, see Table (II)) 


10) D = 22.0 (ft) 


11) Ny 


105 (rpm) 


12) P 30000.0 (hp) 


D 

Using an optimization scheme incorporated in his MVAPDP 
computer program, Vassilopoulos maximized the open water 
efficiency (n,) and designed a propeller with the following 
characteristics: 


1) J = .852 
2) K,, = .242 
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3) 0478 


te 
I 


4) n. = .691 
5) Ap/Ag = .767 


6) bldwt 


7617.2 (lbf) 

By utilizing both the lifting line and lifting surface 
methods in his design procedure, Vassilopoulos' MVAPDP program 
evolved a “constant stress" propeller blade. Consequently, 
the values for (t*/c) and P/D varied non-linearly along the 
propeller radius (R). According to Vassilopoulos, this 
resulted in a minimum weight propeller. The values for P/D 
and (t*/c) are listed in Tables (8) and (10) of his paper. 
From these values, (t*/c) o5p is approximately .040. 

While the propeller represented by Vassilopoulos' design 

? is different, in many aspects (rake, skew, blade section 
aerfoil shape, etc.), from the Wageningen B-Screw Series 
propeller, it does represent a minimum weight propeller that 
has been “matched" to specific design values. Appropriate 


results are summarized in Table (VI). 


D. SOLUTIONS BY COPES/CONMIN 
The propeller selection problem, as stated by equation 
(9.8), is now solved by COPES/CONMIN. One solution variation 
is considered. The following parameters are used: 
1) Temp = 59 (°F) 
1.9905 (1bf-sec*/£t*) 


2) p= 

-5 2 
3) v = 1.2817 x10 (£t"/sec) 
4) -247 (psia) 


Pwatvap = 


5) Pia 14.7 (psia) 

6) wt = .22 

7) Np = 1.025 

8) noscrw = l 

9) Nee = 19.0 (ft) 
10) Z= 6 
ll) promat = 5 (stainless steel, see Table (II)) 
12) D = 22.0 (ft) 


Two problems are examined. Problem 1 specifies the following 


additional parameters: 
1) td = .1725 
2) PE = 21292.6 (hp) 
3) V = 24.24 (knots) 
4) Q, = 1500606.75 (ft-1bf) 
5) N= 105 (rpm) 


Problem 2 specifies the same parameters as: 


1) td = .171 

| 2) P, = 17630.0 (hp) 
3) Vv = 23.0 (knots) 
4) Q, = 1500606.75 (£t-lb£) 
5) N, = 105 (cpm) 


All of the above are initialized in che input section 
(ICALC = 1) of similar versions of SUBROUTINE ANALIZ. There- 
fore, only one version is included in Appendix I. 

The constraints for G,(X3) and Gy (X35) are evaluated by 


SUBROUTINE BLPOW3 which appears in the execution section of 
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SUBROUTINE ANALIZ. Also, note that SUBROUTINE DICNUA has 
been deleted from the execution section, while SUBROUTINE 
WGTCAL has been added. 
1. Programming Details 
All twelve constraints are evaluated (NCON = 12). 
The upper (xa7PPP*) and lower (ape) limits on the design 


variables P/D, Ap/Ap and (t*/c) sep are set to be: 


4 < P/D < 1.4 
4 < Ap/Ag eee res 
003 <= (t*/c) 75R < 50 


These upper and lower limits are speci fied in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (X3,) 
is also assigned on card image F under the field labeled xX. 
The list of card images in Appendix J lists all of the COPES 
control cards used for both problems. These cards also 
specify the locations of the design variables in the common 
block GLOBCM (see Table (III)) as well as the locations of 
the constraints and their boundaries. Further details on 
the COPES control card requirements and the format of each 
card are contained in Reference [7]. 

2. Results 

The outputs from the optimization/analysis, performed 

by COPES/CONMIN, are listed in Appendix K. Results of both 


problems are tabulated in Table (VI). 
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E. DISCUSSION 

Table (VI) presents the results of problems 1 and 2 
along with relevant information from Vassilopoulos' "design". 
Problem 1 attempted to "match" a Wageningen propeller at the 
design point found by Vassilopoulos. The first COPES/CONMIN 
printout in Appendix K indicates that the “match” was achieved 


at P, equal to 21.168.1 (hp) and Pp equal to 28,150.0 (hp) 


E 
(or, Qo = 1500607 (ft~lbf) and N, = 105 (rpm)). These values 


P 
are judged to be close enough to the "Given" values in Table 
(VI). 

It is apparent that the Wageningen propeller does not 
require all of the 30,000 (hp) of delivered horsepower. The 
propeller characteristics (i.e., J, Kn, Ky and nN) for problem 
1 compare very well to Vassilopoulos' values. The expanded 
area ratios (Ap/AQ) are, also, very Similar. Of course, the 
obvious difference is the blade weight (bldwt). The Wageningen 
propeller blade is over five thousand pounds heavier. Does 
this make sense for a minimum blade weight? 

The answer is yes. 

All one has to do is consider the values of (t*/c) for problem 
l and Vassilopoulos' design. Vassilopoulos' "constant stress" 
blade was designed to "absorb" stress up to the allowable 
design limit of 5,400 (psi) (for stainless steel) all along 
the entire propeller radius (R). Table (12) in Reference 

[18] gives further details. The Wageningen propeller blade, 
however, represents an “older” type of blade which was de- 


signed with a linear blade section maximum thickness (t*) 
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distribution. Consequently, it was “overdesigned” for strength 
beyond the 3/10~--4/10 radius (i.e., .3R--.4R) and contains 
excess material. A heavier blade, therefore, results. Note, 
also, that the optimizer did not drive the value of (t*/c) o5p 
to the minimum acceptable value, (t*/C) 3565p min: 

The results of problem 2 show the effect on blade weight 
(bldwt) for a Wageningen propeller when the hull's powering 
requirements (i.e., Pp at V) have been reduced. The weight 
reduction of 2000 pounds is significant. The complete re- 
sults are listed in the second COPES/CONMIN printout in 


Appendix K. 
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TABLE VI 


Design Case No. 3-~Results 


GROUP ITEM VASSILOPOULOS PROBLEM 
l 2 
. PE 21292.6 21292.6 17630 .0 
Given ¥ 24.24 24.24 23.0 
Q 1500607 1500607 1500607 
Np 105 105 105 
Design 
Variable D 22.0 22.0 22.0 
Specified 
P/D * 1.1813 1.0906 
eee AYAy 767 7944 .7742 
(t4/c) ae, .040 .0794 .0681 
Minimize  bldwt aos 12842.6 10464 .7 
Maximize No 691 —_— —_— 
Di im 7 aa 7A 
Restrice RR; ci = 8515 .7722 
(8970) genas = 0691 0681 
J 852 .8290 . 7866 
Kp 242 £2349 .2063 
Other Ky .0478 0448 .0372 
No oe .6915 6950 
blavt 7617.2 sis a 


* 
P/D varies with R 
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X. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The general purpose non-linear optimizer/synthesizer 
COPES/CONMIN has been successfully applied to three typical 
preliminary ship design propeller selection problems in which 
the Wageningen B-Screw Series is used. The formulation and 
programming of each required analysis code (i.e., SUBROUTINE 
ANALIZ) have been made as general as possible to allow the 
designer a broad variety of solution options for solving 
propeller selection problems which can be classified under 
any of the three Design Cases that were considered. The 
analysis codes have been "“modularized" to the extent that 
methodical series data from other propeller series, which 
are available in the polynomial expression format of the 
B-Screw Series, can be easily adapted for powering analysis 
utilizing design optimization methods. 

Further flexibility in the solution to the propeller 
selection problem has been achieved by using COPES/CONMIN 
as the optimizer/synthesizer. The designer has now been 
afforded the additional capability of specifying the design 
variables, the objective functions and the constraints of his 
choice. By solving propeller selection problems in the way 
presented in this thesis, repetitive problem formulation and 
coding have been eliminated. 


There are other advantages to solving propeller selection 


problems specifically with COPES/CONMIN which have not been 
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directly addressed in this study. As stated in Chapter II, 
COPES/CONMIN is capable of performing optimization analyses, 
sensitivity studies, optimum sensitivity studies and optimi- 
zation using approximation techniques. The designer, there- 
fore, can select and perform any of these options, using the 
same analysis codes which have been presented in this thesis. 
While the utilization of a general purpose non-linear 
optimizer in solving propeller selection problems allows the 
designer greater flexibility in the selection procedure, there 
is one important limitation that should be stressed at this 
point. This concerns the question whether or not the solution 
vector, determined by the optimizer, is a "global" optimum. 
As stated in Chapter II, COPES/CONMIN assures that, if a 
feasible solution vector is found, it is, at least, a “local" 
minimum (or maximum). This implies that, for two different 
initial design vectors which are specified in the COPES Con- 
trol Card deck on card image F, the same optimum solution 
may not be determined by the optimizer. Both solutions would 
correspond to minimums (or maximums) of the objective function 
and are, therefore, correct. But, does one or the other 
correspond to the minimum (or maximum) of the entire vector 
design space, i.e., the “global" optimum? For the moment, 
at least, there is no definitive answer to the question. 
Despite this uncertainty, progress in the field of design 
optimization continues to be made. Current developments 


(Ref. 47] will soon allow the designer to have a choice in 
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selecting a specific optimization algorithm froma "library" 
of proven optimization programs which employ the latest 
state-of-theart numerical techniques. Again, using ane 
analysis code, the designer will be able to generate any 


number of optimized solutions for the problem under study. 


B. RECOMMENDATIONS 

For future consideration, it is recommended that the 
automated design and trade-off capability, provided by a 
general purpose non-linear optimizer/synthesizer such as 
COPES/CONMIN, be applied to the more difficult problem of 
propeller design. 

As pointed out in Chapter I, the use of the Wageningen 
B-Screw Series represents a “selection" procedure rather 
than a "design" process. Today, analytical propeller design 
procedures, utilizing lifting line and lifting surface 
theory, are becoming increasingly popular among propeller 
designers. The propeller design, which results from the 
utilization of these analytical methods, is, unquestionably, 
more efficient than the standard series propeller. However, 
these methods require consideration of many more design varia- 
bles in the design process. This appears to be a natural 
application for the use of a general purpose non-linear 
optimizer/synthesizer. 

Here, an analysis code, much larger than those which 


have been presented in this study, could be developed which 


would incorporate the lifting line/lifting surface theory 


for the determination of the propeller performance character- 
istics, the local cavitation numbers and also the calculation 
of the pressure distributions over the blade. These pres- 
sure distributions would be utilized in the strength analysis 
of the blade. This analysis would utilize the finite element 
technique on an appropriately generated mesh model of the 
blade. Having defined the steps for this design procedure 

in the analysis code, the propeller designer now "couples" 
his analysis to the optimizer/synthesizer for determination 
of the optimum design. A massive amount of computer storage 
would certainly be required, but this concept is feasible 


and, in the author's view, is worthy of future consideration. 


C. A FINAL NOTE 

In conclusion, this thesis has demonstrated, in effect, 
another interesting application of the method of design opti- 
mization. The author, in no way, wishes to leave the reader 
with the impression that the techniques of design optimi- 
zation are the "be all--end all" for engineering analysis. 
Design optimization techniques are useful and powerful tools 
that stand to relieve the engineer of the mundane tasks of 
numerical calculations and subsequent graphic plotting. 
But, they are just tools. In the final "analysis", good 
engineering judgment is paramount in their application and 


use « 
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APPENDIX A 


FORTRAN VARIABLE CROSS REFERENCE LIST 


Pwatvap 


Patm 
promat 


Qs 


* 
Rn 75R 


So 


td 
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Fortran Variable 


AEDVAO 
AEAOMN 
WEIGHT 
C75R 
DIA 
DIALIM 
HCL 

RJ 

KQ 

KT 
NOSCRW 
N 

PE 

PD 
PDIVD 
PWATVA 
PATM 
PROMAT 
Qs 
R75R 
sc 


TD 


Symbol 
Temp 

(t*/c) 75R 
Vv (ft/sec) 
Vv (knots) 


wt 


APPENDIX A (CONT.) 


Fortran Variable 
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APPENDIX B 


SUBROUTINE LISTINGS 
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APPENDIX D 


CONTROL CARD IMAGES--DESIGN CASE NO. 1 
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COPES OUTPUT--DESIGN CASE NO. 
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